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EPITAXIAL  SUBLIMATION  METHODS  FOR  THE  STUDY  OF  PSEUDO¬ 
BINARY  SEMICONDUCTOR  ALLOYS 


INTRODUCTION 

This  annual  report  describes  work  performed  under 
ARPA  Order  1597  at  the  Naval  Ordnance  Laboratory.  This  portion 
of  the  ARPA  program  centers  on  surface  properties  affecting 
the  electro-optical  behavior  of  semiconductors  particularly 
the  IV- VI  semiconductors  and  alloys  used  in  infrared 
emitters  and  detectors.  These  surface  properties  are 
dominated  by  chemisorbed  impurities.  (Related  effects  produced 
by  electromigration  of  impurity  ions  dominate  the  electrical 
properties  of  the  semiconductor-oxide  interface.)  In  IV-VI 
device  technology  chemisorption  effects  are  both  strong  and 
diverse,  determining  carrier  populations  in  surface  regions 
and  throughout  film  structures,  producing  sensitization  and  aging, 
and,  as  we  have  now  shown,  affecting  the  quality  of  epitaxial 
growth . 

The  primary  purpose  of  the  NOL  ultra-high  vacuum  experimental 

effort  is  to  clarify  the  nature  of  the  interaction  of  gases, 

primarily  0p  and  H,  with  epitaxial  IV-VI  films.  Previous 
1-6 

studies  of  the  effects  of  chemisorption  on  the  electrical 
properties  of  IV-VI  films  have  been  carried  out  in  relatively 
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modest  vacua  and/or  on  films  that  had  been  exposed  to 

the  atmosphere.  As  a  consequence,  the  role  of  contamination 

in  these  experiments  is  unknown,  and  the  results  ambiguous. 

In  the  NOL  program,  adsorption  studies  are  made  on  films 
which  have  been  grown  and  maintained  under  ultra-high  vacuum 
conditions.  The  Hall  and  resistivity  measurements  are  made 
iH  situ >  and  ambiguities  due  to  contamination  eliminated.  Thus 
far,  we  have  obtained  a  number  of  significant  new  results 
bearing  on  a)  film  growth;  b)  hydrogen  adsorption  on  PbS;  and 
c)  the  oxidation  of  PbS.  In  addition,  we  have  formulated  a  new 
method  of  analyzing  transport  data  which  substantially  reduces 
the  difficulty  of  interpreting  the  behavior  of  two-carrier 
films.  Details  of  this  work  and  its  implication  for  IV-VI 
thin  film  device  technology  are  given  in  sections  1,  2,  3,  and  4. 

The  chemical  state  of  contaminated  IV-VI  surfaces  is 
responsible  for  many  of  the  most  important  electro-optical 
properties,  but  is  not  explained  well  by  bulk  chemistry.  Many 
special  methods  o°  experimental  surface  analysis  have  been 
devised  for  determining  the  chemical  species  present  and  their 
bond  structure.  We  have  built  the  simplest  of  such  surface 
analysis  tools,  an  appearance  potential  spectroscopy 
(APS)  apparatus.  APS  uses  an  electron  beam  to  probe  the  core 
states  of  surface  atoms,  measuring  the  threshold  of  X-ray 
florescence  core  de-excitation.  The  spectrometer  has 
been  assembled  and  the  electronics  is  being  tested. 

Details  are  given  in  Section  5- 
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Our  theoretical  studies  have  made  important  advances  in 
the  understanding  of  the  effect  of  chemisorption  on  the 
electrical  properties  of  semiconductors.  We  have  completed1^ 
a  detailed  calculation  of  carrier  scattering  by  charged 
chemisorbed  ions.  This  takes  into  account  for  the  first  time 
the  actual  three  dimensional  screened  space  charge  potential 
of  the  surface  charge  array,  the  interference  terms  between 
electrons  scattered  from  different  surface  charges,  and  the 
dielectric  image  force;  all  these  effects  are  shown  to  be 
quite  important .  These  results  are  now  being  used  as 
the  basis  of  a  realistic  surface  transport  theory  to  predict 
surface  mobility,  surface  Hall  coefficient,  and  surface 
magnetoresistance.  We  have  also  made  a  breakthrough  in  the 
theory  of  semiconductor  surfaces  with  the  discovery1"7’1 
of  a  new  and  basic  theoretical  relation  between  surface 
scattering  rates  and  chemisorption  thermodynamics.  This  new 
relation  is  based  on  the  structure  factor  of  the  array  of 
chemisorption  charges  on  the  surface.  We  are  following  this 
up  with  calculations  of  the  adsorption  isotherms  of  various 
gases  on  semiconductor  surfaces.  Progress  has  also  been  made 
on  the  quantum  theory  of  the  dispersive  dielectric  image  force, 
which  we  have  shown  to  be  important  in  surface  mobility,  etc., 
and  which  is  well-known  to  strongly  affect  surface 
recombination  rates,  heterojunction  tunneling  etc. 

Details  of  this  work  are  given  in  sections  6,  7,  8,  9,  and  10. 
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RESEARCH  PROGRESS 


1.  Film  Growth  (Experimental:  R.  N.  Lee) 

We  find  that  the  room  temperature  mobilities  In  our  ultra- 
high  vacuum  grown  PbS  films  are  comparable  to  or  better  than 
those  In  the  best  films  grown  In  less  elaborate  vacuum 
systems.  The  measured  room  temperature  mobilities  range  from 
approximately  500  cm2/v  lt-sec.  to  as  high  as  750  cm2/volt-sec . 
The  higher  values  approach  the  highest  mobilities  that  have 
been  observed  in  bulk  material. 

We  have  considered  the  possibility  that  the  high  apparent 
mobilities  are  duo  to  the  formation  of  a  surface  layer  of  lead. 
The  growth  of  such  layers  has  been  reported  for  sublimated 
PbS  films , 7  and  would  lead  to  an  anomalously  high  apparent 
mobility  for  the  PbS.  The  measured  mobility  for  such  a  film 
would  be  the  resultant  of  the  carrier  mobility  in  the  lead 
layer  and  the  mobility  in  the  underlying  PbS  layer.  This 
possibility  may  be  rejected,  however,  on  the  basis  of  the 
adsorption  results  to  be  discussed  later  in  this  report.  The 
adsorption  data  clearly  indicate  that  there  is  only  one  carrier 
type  (i.e.  only  one  mobility)  in  the  film,  and  we  may  have 
confidence  in  the  high  measured  values  of  the  mobilities. 

The  superior  quality  of  the  ultra-high  vacuum  grown  films 
gives  us  significant  new  insight  into  the  mechanisms  of  growth 
and  strain  relief  on  NaCl  substrates.  Epitaxial  IV- VI  films 
are  commonly  grown  in  vacua  of  ^10-5  Torr  with  water  vapor 
comprising  90 %  or  more  of  the  background  gas.  A  popular 
conjecture  has  been  that  the  water  is  necessary  for  good  film 
growth  1 i  that  it  allows  the  NaCl  to  flow  plastically  and/or 

provides  a  buffer  between  the  substrate  and  the  film.  This 
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Idea  must  now  be  discarded,  for  the  water  vapor  in  the  ultra- 
high  vacuum  system  Is  completely  negligible  and  the  substrates 
were  thoroughly  degassed  prior  to  deposition  of  the  films. 

These  results  suggest  that  substantial  improvements  in 
device  quality  could  be  obtained  through  the  use  of  ultra- 
high  vacuum  during  device  fabrication.  We  note,  in  fact,  that 
bakeout  of  the  vacuum  system  was  one  of  the  procedures  employed 

r\ 

by  Holloway  et  al°  in  obtaining  their  very  high  mobility  films, 
and  we  suggest  that  this  procedure  is  of  greater  importance 
than  has  previously  been  appreciated. 

2.  Adsorption  of  Atomic  Hydrogen  (Experimental:  R.  N.  Lee) 

One  of  the  major  results  of  recent  work  by  Zemel  and  co- 
workers^  >5  at  the  University  of  Pennsylvania  has  been  that 
the  electrical  properties  of  epitaxial  lead  chalcogenide 
films  are  highly  sensitive  to  atomic  hydrogen.  The  hydrogen 
effects  observed  were  completely  reversible,  and  Zemel  has 
suggested^  that  hydrogen  adsorption  may  provide  a  non¬ 
destructive  measure  of  the  lead  vacancy  concentration  in  lead 
salt  films.  This  exciting  possibility  presupposes  a  model 
in  which  the  hydrogen  occupies  lead  vacancy  sites.  5  The 
circumstances  of  the  previous  studies  were  such,  however, 
that  the  hydrogen  effect  could  be  associated  with  the  presence 
of  oxygen  on  the  film  rather  than  with  the  lead  vacancies,  and 
it  is  of  some  importance  to  examine  the  interaction  of  hydrogen 
with  the  lead  salts  in  the  absence  of  cxygen  contaminatior . 

For  this  reason,  we  have  studied  the  adsorption  o:  hydrogen 
on  ultra-high  vacuum  grown  PbS  films.  The  adsorption  effect? 
were  observed  by  means  of  in.  situ  measurements  of  the  Hall 
coefficient  and  the  resistivity. 
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The  ultra-high  vacuum  studies  have  verified  the  major 
features  of  hydrogen  adsorption  which  were  previously  reported. 

That  is,  molecular  hydrogen  has  no  observable  effect  on  the 
electrical  properties;  whereas  atomic  hydrogen  adsorbs 
reversibly  to  produce  relatively  large  changes  in  the  carrier 
population.  Atomic  hydrogen  acts  as  a  donor,  and  it  is 
incorporated  into  the  film  without  producing  significant 
changes  in  the  electron  mobility. 

Figure  1  is  an  RH  versus  p  plot  which  resulted  from 
exposure  to  atomic  hydrogen.  This  data  is  somewhat  complicated, 
and  the  lettered  key  points  in  the  plot  are  identified  in 
Table  1.  The  major  points  to  be  derived  from  this  plot  are 
that  the  hydrogen  increased  the  electron  population  by  up  to 
^•5  3n  cm  and  that  the  data  points  fell  along  a  constant 
mob ill tv  line  as  long  as  the  temperature  variation  remained 
negligible . 

In  our  investigations,  the  atomic  hydrogen  was  generated 
by  a  beam  of  70  eV  electrons  in  the  ionization  stage  of  the 
residual  gas  analyzer  mounted  on  the  experimental  chamber, 
and  the  atomic  hydrogen  so  produced  had  to  make  several 

collisions  with  the  chamber  walls  before  reaching  the  sample.  « 

The  effective  partial  pressure  of  atomic  hydrogen  at  the 

sample  is  thus  very  difficult  to  estimate.  It  k,s  assumed, 

however,  that  the  amount  of  atomic  hydrogen  at  the  sample 

was  proportional  to  the  total  hydrogen  pressure  in  the  system 

as  measured  by  a  Bayard-Alpert  gauge.  It  should  be  pointed 

out  that  the  metal  walls  in  the  baked  ultra— high  vacuum 

system  pump  atomic  hydrogen  much  more  efficiently  than  the 
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TABLE  I 


KEY  TO  FIGURE  1 


Key  Letter 

Event 

a 

The  residual  gas  analyzer  was  turned  on,  result¬ 
ing  in  an  increased  pressure  and  production  of 
atomic  hydrogen  at  a  pressure  estimated  to  be 
no  higher  than  5  x  10“10  Torr. 

b 

After  six  hours,  the  RGA  was  turned  off. 

c 

Hydrogen  gas  was  admitted  to  the  experimental 
chamber  at  a  pressure  of  5  x  10~3  Torr. 

d 

A  baffled  ionization  gauge  was  turned  on, 
producing  a  small  amount  of  atomic  hydrogen 
at  the  sample. 

e 

The  RGA  ionizer  was  turned  on  to  generate 
atomic  hydrogen  more  efficiently. 

f 

This  extreme  data  point  was  taken  immediately 
after  turning  on  the  RGA.  Data  points  taken 
during  the  subsequent  decrease  in  pressure 
fell  along  f-g. 

g 

The  total  hydrogen  pressure  was  stabilized 
at  5  x  10“5  Torr. 

h 

The  pressure  was  again  allowed  to  fall,  and 
data  points  along  h-k  were  taken. 

k 

The  pressure  was  again  raised  to  the  10“^  Torr 
range.  Subsequent  variations  in  pressure  and 
temperature  resulted  in  the  data  points  lying 
between  e-g  and  h-k. 
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unbcikca  glass  walls  of  the  chambers  used  In  the  previous 
studies,  and  the  effective  pressure  of  atomic  hydrogen  was 
undoubtedly  much  lower.  This  difference  In  pressure  regimes 
Is  reflected  In  the  data. 

The  lower  atomic  hydrogen  pressures  used  In  the  present 
study  resulted  in  the  observation  of  some  previously  unreported 
aspects  of  the  adsorption.  We  have  found  that  the  hydrogen 
uptake  at  these  pressures  Is  a  function  of  pressure  over  a 
pressure  range  of  more  than  four  orders  of  magnitude.  That 
Is,  when  a  constant  hydrogen  pressure  was  established,  the 
electron  population  quickly  stabilized  at  a  steady-state 
value  which  depe*  "'d  on  the  pressure.  A  plot  of  the  steady 
state  electron  population  versus  total  hydrogen  pressure  Is 
presented  In  Fig.  2. 

The  kinetics  of  the  hydrogen  adsorption  could  not  readily 
be  determined  in  the  ultra-high  vacuum  studies  due  to  the  well- 
known  'pump-down'  effect.  As  soon  as  the  electron  beam  was 
turned  on,  the  pumping  of  the  atomic  hydrogen  produced  caused 
a  rapid  drop  in  the  total  pressure,  and  it  was  usually  many 
minutes  before  a  steady  state  was  achieved.  During  this  time, 
the  effective  pressure  of  atomic  hydrogen  at  the  sample  was 
not  proportional  to  the  total  hydrogen  pressure  and  it 
certainly  was  not  constant. 

When  the  hydrogen  pressure  was  suddenly  changed  from  one 
constant  value  to  another,  the  new  steady-state  value  of  the 
electron  population  was  usually  established  more  quickly  than 
child  be  followed  with  our  point-by-point  measurement  technique. 


9 


LOG  P(H2)  (TORR) 


We  are  now,  changing  our  instrumentation  to  provide  a  continuous 
recording  of  both  the  resistivity  and  Hall  voltages,  however, 
and  the  kinetics  of  these  transitions  should  be  available  in 
the  future . 

The  adsorption  data  have  all  been  plagued  with  the  effects 
of  temperature  drifts.  Although  the  effects  are  quite  large, 
no  mention  of  these  difficulties  has  been  made  in  recent 
reports  of  similar  transport  studies.  The  problem  arises  as 
a  consequence  of  the  T~5/2  temperature  dependence  of  the 
mobility  in  the  lead  salts.  This  temperature  dependence 
implies  that  air10  temperature  change  at  room  temperature 
produces  a  0.8%  change  in  the  resistivity,  and  this  result 
is  entirely  consistent  with  our  observations.  Normal  short 
term  temperature  variations  in  our  laboratory  are  of  the  order 
of  a  degree,  and  a  five  degree  drift  over  a  four 
hour  period  is  not  unusual.  These  temperature  fluctuations 
correlate  well  with  the  scatter  of  the  resistivity  data,  and 
are  primarily  responsible  for  the  relatively  large  scatter  in 
the  data  of  Fig.  1.  The  horizontal  bar  in  the  lower  part  of 
this  figure  indicates  the  resistivity  change  predicted  by  the 
T  5/2  dependence  for  a  two  degree  temperature  drift.  The 
oxygen  adsorption  data  which  are  discussed  later  in  this 
report  were  obtained  during  periods  when  the  ambient  temperature 
was  relatively  stable,  and  the  scatter  of  these  data  is 
significantly  smaller. 
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3*  Oxygen  Adsorption  (Experimental:  R.  N.  Lee) 

Our  studies  of  the  chemisorption  of  oxygen  on  PbS 
constitute  the  first  examination  of  the  oxidation  of  a 
lead  chalcogenide  to  be  carried  out  by  means  of  transport 
measurements  in  the  very  low  pressure  range.  We  have  made 
controlled  oxygen  exposures  at  pressures  ranging  from  10"® 
Torr  up  to  10“2  Torr  with  partial  pressures  of  background 
gases  less  than  10-10  Torr.  Previous  studies1-6  were 
carried  out  with  background  pressures  of  10-^  Torr  or  higher 
and  with  oxygen  exposures  primarily  provided  by  exposing 
the  films  to  atmosphere.  Our  results  include  a  number  of 
previously  unobserved  features  of  the  oxidation  of  lead 
chalcogenides  and  show  that  some  apparent  discrepancies 
between  previous  studies  are  the  result  of  differing  vacuum 
conditions . 

The  effects  of  exposing  PbS  to  oxygen  at  pressures  up  to 

10-2  Torr  are  shown  in  the  Hall  coefficient  versus  resistivity 

plot  of  Pig.  3.  Oxygen  acts  as  an  acceptor  in  PbS  (as  it  does 

in  all  of  the  lead  salts)  and  exposures  over  this  pressure 

range  produced  a  maximum  change  in  the  electron  population  of 

2.6  x  lO1^  cm-2.  This  corresponds  to  a  23#  reduction  in  the 

number  of  electrons  in  the  conduction  band.  Although  the  RTJ 

n 

versus  p  plot  does  not  reveal  it,  part  of  this  oxygen  uptake 
was  reversible  and  pumping  to  low  pressures  always  provided  a 
partial  restoration  of  the  electron  population. 

We  note  in  particular  that  the  data  points  of  Fig.  3  lie 
very  precisely  along  a  straight  line  whose  slope  gives  a 
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mobility  of  63O  cm^/volt-sec .  This  degree  of  linearity  could 
not  be  achieved  if  the  film  contained  more  than  one  type  of 
carrier,  for  that  would  require  unrealistic  variations  in  the 
carrier  populations  and/or  mobilities.  The  data  of  Pig.  3 
therefore  constitute  firm  evidence  that  the  film  cor. '-ains  only 
a  single  carrier  type  having  a  constant  mobility. 

The  very  general  results  just  cited  are  quite  similar  to 
the  results  previously  reported  for  the  lead  chalcogenides , 1~r 
and  our  observations  remove  a  number  of  the  fundamental 
ambiguities  implicit  in  the  previous  studies.  That  is,  we  have 
for  the  first  time  eliminated  the  possibility  that  surface 
oxides  or  high  background  pressures  of  active  gasses  are 
responsible  for  major  features  of  the  adsorption  such  as  the 
large  changes  in  carrier  population  or  the  reversibility  of 
the  adsorption.  However,  we  have  also  obtained  a  number  of 
new  results  which  indicate  that  the  oxidation  of  the  lead 
salts  may  be  appreciably  more  complicated  than  previous 
transport  studies  would  suggest. 

Our  observations  show  that  the  oxidation  of  PbS  comprises 
at  least  four  processes  which  can  be  distinguished  on  the 
basis  of  rate  or  reversibility.  Our  observations  have 
revealed  two  reversible  processes.  One  of  these  takes  place 
at  a  very  rapid  rate  and  saturates  in  less  than  two  minutes 
at  pressures  as  low  as  2  x  10“h  Torr.  The  second  of  these 
processes  proceeds  at  a  slower  rate  and  saturates  at  an 
oxygen  uptake  which  is  dependent  upon  the  oxygen  pressure. 

We  have  also  observed  that  there  are  one  or  more  processes 
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which  result  in  an  irreversible  oxygen  uptake.  The  rate 
(or  rates)  of  the  irreversible  adsorption  have  not  been 
directly  observable,  but  it  has  been  found  that  the 
Irreversible  uptake  saturates  in  an  hour  or  less  at  a 
pressure  of  6  x  1C  ^  Torr.  We  deduce  from  results  in  the 
literature  >10>H  that  there  must  be  at  least  one  more 
extremely  slow  process  which  produces  the  surface  oxide 
layer  which  is  known  to  form  when  PbS  is  exposed  to  the 
atmosphere . ^ u > H 

Each  of  the  three  chemisorption  processes  we  observe  in 
the  low  pressure  regime  can  be  associated  with  a  specific 
component  of  the  total  change  in  the  electron  population. 

The  Irreversibly  adsorbed  oxygen  is  always  responsible  for  a 
J-arge  fraction  of  the  total  electron  population  change  and, 
of  course,  once  saturation  is  reached  this  change  is  permanent. 
In  the  case  of  the  reversibly  adsorbed  oxygen,  the  rapid 
process  is  responsible  fo~  only  a  small  fraction  of  the  total 
population  change,  while  the  slower  process  accounts  for  the 
balance.  In  a  typical  exposure  run  at  10  2  Terr,  irreversible 
adsorption  removes  lO1^  electrons  per  cm2,  the  rapid  reversible 
process  removes  2  x  1012  electrons  per  cm2  and  the  slower 
reversible  process  removes  1.4  x  lO1^  electrons  per  cm2  for 
a  total  electron  population  decrease  of  2.6  x  lO1^  cm"2. 

These  results  are  illustrated  by  Figs.  4,  5  and  6.  The 
curves  of  Fig.  4  were  both  obtained  under  the  initial 
condition  that  the  film  contained  no  reversibly  adsorbed 
oxygen,  and  the  initial  rapid  uptake  is  readily  apparent.  For 


15 


TIME  (MINUTES) 


I 


TIME  ■  MINUTES) 


curve  (a),  the  film  had  received  no  previous  oxygen  exposure 
and  consequently  contained  no  adsorbed  oxygen  at  all.  Curve 
(b)  was  obtained  from  the  same  film  after  the  irreversible 
adsorption  had  been  saturated  and  all  of  the  reversibly 
adsorbed  oxygen  pumped  off.  We  note  that  the  amount  of 
rapidly  adsorbed  oxygen  was  the  same  in  both  cases  even 
though  the  oxygen  pressure  was  2  x  10“ft  Torr  in  one  instance 
and  1.5  x  10“^  Torr  in  the  other.  We  observe,  as  well,  that 
the  subsequent  slow  uptake  occurred  at  a  faster  rate  in  the 
initial  exposure  in  spite  of  the  fact  that  the  pressure  was 
higher  than  in  the  later  exposure.  Thus,  the  slow  uptake 
rate  in  curve  (a)  is  a  composite  of  the  rates  of  the 
irreversible  and  the  slow  reversible  processes. 

The  curve  of  Fig.  5  is  a  continuation  of  Fig.  4(b).  It 
illustrates  that  the  slow  process  saturates  at  a  constant 
oxygen  pressure  and  a  steady  state  electron  population  is 
achieved.  The  slow  process  clearly  occurs  at  a  pressure 
dependent  rate,  and  we  note  that  saturation  was  achieved 
relatively  rapidly  alter  the  pressure  was  increased  to 
2.4  x  10“^  Torr. 

In  Fig.  6,  we  see  an  example  of  the  desorption  that 
occurs  when  the  oxygen  is  pumped  out  of  xperimental 

chamber.  In  this  case,  oxygen  was  re-a  i‘  i  before  all  of 
the  reversibly  adsorbed  oxygen  was  removed,  and  an  initial 
rapid  uptake  was  not  observed.  Evidently,  the  rapid  process 
remains  saturated  until  all  of  the  oxygen  associated  with  the 
slow  process  has  been  removed.  The  successive  pressure 
increases  indicated  in  Fig.  6  show  that  the  saturation  of  the 
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slew  process  is  pressure  dependent. 

These  results  were  obtained  at  oxygen  pressures  which 
were  five  to  ten  orders  of  magnitude  lower  than  those  used 
in  previous  transport  studies  of  oxygen  chemisorption  on 
lead  chalcogenides , and  a  direct  comparison  with  the 
high  pressure  results  cannot  validly  be  made.  It  is 
unlikely,  for  example,  that  what  we  term  the  rapid  process 
has  ever  been  seen  before.  We  believe,  in  fact,  that  the 
fast  initial  rate  reported  by  Egerton  and  Juhasz^  is  actually 
a  high  pressure  version  of  the  slow  rate  reported  here,  and 
that  their  slow  process  is  related  to  the  growth  of  a  surface 
oxide . 

Our  results,  when  viewed  in  the  light  of  the  known 
properties  of  PbS,  raise  serious  questions  about  the  validity 
of  the  diffusion  models  which  have  been  proposed  for  the 
chemisorption  of  oxygen  on  PbS.2»^5  if  we  assume  that  the 
irreversible  adsorption  is  associated  with  surface  defect 
structure,  then  a  diffusion  model  would  require  that  the 
rapid  uptake  correspond  to  ionic  adsorption  at  the  surface 
and  the  slow  uptake  correspond  to  diffusion  into  the  bulk. 
However,  this  would  seem  to  leave  no  mechanism  to  explain  the 
very  slow  uptake  that  has  been  observed  to  affect  the  carrier 
concentration  in  previous  transport  studies.  Furthermore, 
it  would  seem  much  more  likely  that  a  surface  mechanism 
would  display  the  observed  reversibility  and  pressure 
dependence . 

It  would  seem  that  any  proposed  model  would  have  to  be 
consistent  with  the  formation  of  a  surface  oxide  in  order  to 
be  entirely  plausible,  and  the  diffusion  model  does  not 


20 


appear  to  be  consistent  with  the  known  facts.  Low-energy 
electron  diffraction  (LEED)  studies^  tell  us  that  a  48 
Torr-hour  exposure  to  oxygen  results  in  approximately  0.1 
monolayer  of  amorphous  oxide  on  the  surface  and  that  exposure 
to  atmosphere  produces  a  thick  surface  layer  of  amorphous 
oxide.  In  previous  transport  studies,  a  48  Torr-hour  exposure 
was  exceeded  in  the  first  four  minutes  at  atmospheric  pressure, 
and  the  entire  oxygen  uptake  observed  in  this  time  would  have 
to  remain  on  the  surface  in  order  to  be  consistent  with  the 
LEED  results.  Also,  the  diffusion  model  makes  no  provision 
for  the  formation  of  a  heavy  surface  oxide.  The  oxide  formation 
would  necessarily  affect  the  carrier  population  in  the  film 
and  it  would  undoubtedly  contribute  to  the  slow  change  in  carrier 
population  that  has  been  reported  to  occur  in  air. 

The  difficulties  in  explaining  oxygen  adsorption  on  the 
lead  salts  have  arisen  primarily  because  of  the  very  large 
charge  exchange  which  has  been  associated  with  the  observed 
changes  in  carrier  population.  The  primary  virtue  of  the 
diffusion  model  is  that,  by  distributing  this  excess  charge 
throughout  the  bulk,  it  avoids  the  necessity  of  extreme 
band-bending.  We  believe,  however,  that  surface  processes 
exist  which  can  explain  all  of  the  data  at  hand  without 
creating  the  extreme  band-bending  which  is  so  objectionable. 

We  are  currently  investigating  a  category  of  surface  processes 
which  may  have  the  desired  characteristics. 
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4.  Analysis  of  Transport  In  Two-Carrier  Films  (Experimental, 


Theoretical:  R.  N.  Lee) 

One  of  the  prominent  results  of  chemisorption  studies  on 
IV-VI  films  has  been  that  the  carrier  mobilities  remain 
unchanged  by  relatively  large  variations  in  the  carrier 
populations.1-6  Interpretation  of  transport  data  is 
considerably  simplified  when  the  carrier  mobilities  are  constant 
and  this  feature  of  the  data  has  been  explicitly  used  in 
previous  interpretations  of  the  effects  of  chemisorption  on 
the  electrical  properties  of  IV-VI  films."'"  6  Previous  analyses 
have  all  assumed  specific  carrier  concentration  profiles 
in  the  films,1-6  however,  and  the  fundamental  simplifications 
provided  by  the  constant  carrier  mobilities  have  not  been  used 
to  the  best  advantage. 

We  will  now  show  that  when  the  carrier  mobilities  are 
constant,  the  transport  data  provide  a  clear  physical  picture 
of  the  changes  occurring  in  a  two-carrier  film.  We  will 
show  that  transport  data  can  tell  us  the  number  of  carriers  of 
each  type  in  the  film  and  give  us  the  change  in  the  number 
of  ionized  dopants  as  well  as  the  degree  of  spatial  overlap 
of  the- two  carrier  types.  This  is  the  maximum  information  which 
may  be  derived  from  the  transport  data,  and  it  is  obtained 
without  introducing  assumptions  as  to  the  carrier  concentration 
profiles  within  the  film. 

In  general,  the  conductivity  a  and  the  Hall  coefficient 
in  two-carrier  layered  films  are  governed  by  the  relatively 
complicated  Petritz12  equations  (we  take  the  carriers  to  be 
electrons  and  holes  for  purposes  of  discussion). 
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ad  =  e/d  y  (z)p(z)dz  +  e/d  y  (z)n(z)dz  , 
o  p  o  n 


(1) 


RHa2d  =  e/d  y  2(z)p(z)dz  -  e/d  y  2(z)n(z)dz  , 
o  p  o  n 

where  d  is  the  film  thickness  and  z  is  the  distance  from  the 
surface . 

When  the  carrier  mobilities  are  constant,  these  relations 
simplify  to 


where 


od  =  eupCp  + 

ey  C  , 

n  n  5 

(3) 

RHo2d  =  eup2Cp  - 

p 

■  ey  C  , 
pn  n  5 

(4) 

Cp  =  /d  p  (z  )dz  and  Cn  =  Jd  n(z)dz 


are  the  hole  and  electron  populations  (total  number  in  the 
film  per  unit  film  area),  respectively. 

We  stress  that  the  use  of  the  carrier  populations  rather 
than  the  carrier  concentrations  avoids  the  problems  associated 
with  guessing  the  carrier  concentration  profiles  p(z)  and 
n(z)  in  (1)  and  (2).  In  fact.  It  is  clear  from  the  equations 
that  Rh  and  a  provide  no  information  whatever  about  the  spatial 
distribution  of  the  carriers  within  the  film. 

Inversion  of  (3)  and  (4)  gives  the  carrier  populations 
directly  in  terms  of  the  experimentally  determined  quantities. 
That  is. 
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(6) 


(RH°%> 


U  +U 
p  Mn' 


r  =  £d_  (VRHg) 
n  ep  (p  +y  ) 
n  VMp  Mn' 


(7) 


The  maximum  possible  information  is  extracted  from  the 
resistivity  and  Hall  data  when  we  relate  the  carrier  populations 
to  the  number  of  ionised  dopant  atoms  added  to  the  film. 

This  is  accomplished  by  invoking  electrical  neutrality  of  the 
film. 


q  _N_  +  eC 
D  D  P 


(8) 


where  C°  and  C°  are  the  initial  carrier  populations  and 
is  the  charge  on  the  ionized  dopant  atom.  We  may  then  use 

6,  7  and  8  to  write  ND  in  terms  of  experimentally  determined 
quantities ; 


nd  - 


— - - [(uSa2 

%unvp  H  0 


-  RHo  )  +  (vn  - 


y(°0  •  o)j 


where  R°  and  oq  are  the  values  obtained  for  the  initial 
condition  of  the  film. 

Thus,  we  see  that  measurements  of  R  and  a  provide  a 

direct  physical  picture  of  the  film.  Plots  of  C  and  C 

p  n 

versus  or  of  C  ,  and  versus  other  experimental 
variables  (e.g.  time  or  gas  pressure)  present  the  measurements 
in  a  highly  intuitive  way.  Such  population  plots  offer  great 
advantages  over  the  commonly  used  R^  versus  p  plots'1-  5"  in 
the  interpretation  of  the  data,  for  the  relationship  between  the 

form  of  the  R^  versus  p  curve  and  the  behavior  of  the  carriers 
is  not  usually  readily  apparent  in  a  two-carrier  regime. 
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The  transverse  magnetoresistances  may  also  be  expressed 
in  terms  which  simplify  interpretation.  It  has  been  observed 
that  there  is  an  anisotropy  in  the  transverse  magnetoresistance 
in  IV-VI  films  that  is  a  function  of  gas  adsorption," 
but  the  basic  reasons  for  this  anisotropy  have  not  been 
adequately  explaineo.  in  the  context  of  the  observed  constant  carrier 
mobilities . 

The  expression  obtained  by  Petritz12  for  the  transverse 
magnetoresistance  with  the  magnetic  field  perpendicular 
to  the  film  is 


(10) 

Although  Petritz  did  not  derive  an  expression  for  the 
case  of  the  magnetic  field  parallel  to  the  film,  McLane  has 
recently  obtained  an  expression  for  this  quantity.^ 
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If  we  again  assume  a  two-carrier  system  with  constant 
mobilities,  equations  10  and  11  reduce  to 


£e_ 

pH 


y  u  e 

=  _£  n 


(ad) 


2(yp 


+  y  )2C  C 
n  n  p 


(12) 


and 


25 


(13) 


Written  In  this  manner,  the  equations  Immediately 
clarify  the  physical  reasons  for  the  anisotropy  in  the 
transverse  magnetoresistance.  The  perpendicular  field 
magnetoresistance  depends  only  on  the  total  number  of  carriers 
of  each  type  in  the  film,  and  is  independent  of  the  carrier 
concentration  profiles.  The  parallel  field  magnetoresistance 
is,  however,  a  measure  of  the  spatial  overlap  of  the  two 
carrier  types.  The  parallel  magnetoresistance  is  therefore 
expected  to  be  less  than  the  perpendicular  magnetoresistance 
in  layered  films  on  very  general  grounds. 

It  is  important  to  note  that  these  expressions  are 
both  zero  if  there  is  only  one  carrier  type  and  tnat  this  is 
an  artifact  of  the  constant  mobility  assumption.  In 
actuality,  of  course,  there  will  be  a  magnetoresistance  in 
a  homogeneous  single-carrier  film  due  to  the  energy  dependence 
of  the  scattering  time.  Such  an  energy  dependence  is  inconsistent 
with  the  strict  insistence  that  the  carrier  mobilities  be 
constant,  however,  for  it  would  imply  that  a  variation  in  carrier 
concentration  would  result  in  a  change  in  mobility. 

If  we  had  maintained  a  distinction  between  the  conductivity 
mobility  <p>,  the  Hall  mobility  <y2>1/,c,  and  the  magnetoresistance 
mobility  Eq.  12  would  also  contain  terms  in  C^2,  Cj2, 

and  which  go  to  zero  when  the  three  mobilities  are  equal. 

(i.e.  when  there  is  no  energy  dependence  of  the  scattering  time.) 
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Analogous  energy  dependent  terms  would  appear  In  Eq.  (13). 

In  practice,  these  terms  will  not  be  zero  and  they  cannot 
easily  be  subtracted  from  the  data.  Thus,  equation  12  may 
not  be  used  quantitatively  with  any  confidence. 

The  size  of  the  energy  dependent  coefficients  in  the 
magnetoresistances  may  be  estimated  by  observing  the  variation 
in  the  measured  magnetoresistances  in  the  single  carrier  regime. 
If  the  magnitudes  and  dependence  on  carrier  popular  ion  in 
this  region  are  small  compared  with  the  variation  in  the 
two  carrier  regime,  it  is  probably  safe  to  make  a  qualitative 
interpretation  of  the  behavior  in  the  latter  region. 

The  population  plot  of  Fig.  7  provides  an  illustration  of 
the  great  utility  of  such  plots.  In  this  hypothetical 
example,  we  presume  that  an  n-p  junction  device  is  fabricated 
by  ion  implanting  acceptors  in  an  n-type  film  and  that 
the  device  fails  after  a  certain  period  of  operation.  We 
see  that  the  population  plot  clearly  identifies  the  failure 
mode  as  being  due  to  electromigration.  The  line  segments 
AB  and  A'B'  correspond  to  the  formation  of  the  junction. 

As  the  donors  are  implanted,  the  electron  population 
decreases  and  the  hole  population  increases  until  imj lantation 
is  terminated  (points  B  and  B').  Operation  of  the  device 
results  in  the  reduction  of  both  populations  by  equal  amounts, 
but  the  amount  of  ionized  dopant  remains  unchanged  during 
this  time  (line  segments  BC  and  B'C').  Such  behavior  of  the 
carrier  and  dopant  populations  could  only  be  the  result  of  the 
diffusion  of  dopant  across  the  junction  and  the  consequent 
mutual  annihilation  of  electrons  and  holes  as  the  junction 
becomes  smeared  out.  Confirmation  of  this  interpretation 
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would  be  provided  by  an  increase  in  the  parallel  magnetoresistance 
during  the  decay  of  the  carrier  populations. 

Alternate  failure  modes  would  have  produced  their  own 
distinctive  population  plots.  For  example,  failure  due 
to  poisoning  by  donor  impurities  diffusing  from  the  surface 
of  the  device  would  result  in  a  population  plot  of  the  type 
illustrated  in  Fig.  8.  These  two  failure  modes  would 
also  produce  different  versus  p  curves,  it  is  true, 
but  only  the  population  plots  clearly  identify  the  mechanism 
at  work. 

The  many  applications  which  can  be  envisioned  for  this 
analytical  technique  indicate  that  it  will  prove  to  be 
invaluable  both  as  a  tool  for  device  diagnostics  and  as 
an  aid  in  the  investigation  of  the  fundamental  phenomena 
underlying  device  technology.  This  work  is  being 
prepared  for  publication. 

5.  Appearance  Potential  Spectroscopy  of  Surfaces  (Experimental: 

B.  V.  Kessler  and  R.  N.  Lee) 

a.  The  method:  It  is  important  to  know  the  chemical 

constituents  of  the  PbS  surface  and  how  they  are  chemically 

bonded,  answering  such  questions  as:  Is  oxygen  present  and 

what  sort  of  Pb;0;S  structures  are  present?  Is  free  lead 

1 3 

present,  etc?  Appearance  potential  spectroscopy  J  (APS)  is  an 
extremely  simple  method  of  surface  chemical  analysis,  not  requiring 
complicated  electron  energy  analysers.  Only  the  soft  X-ray  de¬ 
excitation  of  surface  atoms,  bombarded  by  an  Incident  electron  beam, 
is  measured.  (More  accurately,  ne  measures  the  slope 
of  the  soft  X-ray  yield  versus  electron  beam  energy,  using 
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modulation  techniques.)  Specific  chemical  constituents 
have,  of  course,  specific  core  levels,  and  the  chemical 
shift  of  these  provides  information  about  the  chemical 
environment  of  these  atoms.  At  threshold  for  excitation 
of  an  inner  core  electron  to  the  Perm!  level,  there  is  a 
small  but  sudden  increase  in  the  soft  X-ray  yield  which 
then  produces  a  large  signal  (A  Yield/A  Energy). 

APS  has  been  attempted111  with  III-V  semiconductors  but 
considered  impractical  there  because  the  APS  signals  obtained 
were  several  orders  of  magnitude  smaller  than  with  transition 
metals.  The  reason  for  this  small  signal  in  semiconductors 
is  that  APS  signals  depend  on  the  square  of  the  state 
density  Np  at  the  Fermi  level.  However,  in  the  degenerate 
IV-VI  semiconductors  Np  is  much  higher  than  in  non-degenerate 
semiconductors  so  that  the  method  may  well  be  feasible. 

APS  will  be  applied  first  to  the  TiNi  alloys,  which  have 
high  Np ,  to  test  the  APS  spectrometer,  which  we  have  built, 
and  to  examine  chemical  shifts  at  the  transition  temperature 
of  this  material,  a  second  order  phase  transjtion  <mich 
is  thought  to  be  non-struc tural . 

Although  intended  primarily  as  a  surface  diagnostic  tool, 
the  APS  technique  is  new  enough  to  be  of  some  interest  in 
itself.  The  reasons,  for  example,  why  APS  signals  are  not 
seen  with  some  4d  and  5d  transition  metals  and  alloys  are 
not  yet  understood.  It  also  remains  to  work  out  the  application 
of  the  APS  method  to  materials  used  for  infrared  device 
technology . 
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b.  The  chemical  shift:  With  APS,  surface  impurities 
can  be  identified  with  good  sensitivity.  APS  also 
measures  the  chemical  shift  from  which  one  can 
deduce  information  about  the  structure  of  the  surface. 

This  is  d  :>ne  by  using  an  inner-core-electron  energy  level 
as  a  ’window'  in  which  to  sense  the  charge  in  the  valence 
electrons  of  a  surface  atom  since  a  change  in  the  valence 
state  of  a  surface  atom  results  in  a  variation  of  the  inner 
core  energy  level  and  is  revealed  in  the  APS  electronic  spectra 
The  origin  of  the  chemical  shift  can  be  illustrated 
in  an  idealized  example.  For  an  atom,  consider  only  the 
deepest  inner  K-shell  electron  and  the  outermost  valence 
electron,  ignoring  the  intermediate  electrons.  For  simplicity 
assume  the  valence  electron  has  a  spherical  charge  distribution 
Figure  1  then  delineates  the  electrostatic  potentials  which 
result.  Note  that  the  spatial  extent  of  the  valence  electron 
is  of  the  order  of  1  8  whereas  the  inner  core  electron 
radius  is  much  smaller,  e.g.  can  be  as  small  as  0.01  8  for  a 
high  Z  ion.  Now  it  is  the  valence  electronic  charge 
distribution  which  partakes  in  chemical  interactions  with 
neighboring  ions.  The  effect  therefore  of  a  variation  in  say 
the  magnitude  of  the  assumed  spherical  valence  charge  is  to 
translate  the  nuclear  attractive  potential  which  binds  tne 
inner  core  electron  upward  or  downward  as  in  Fig.  9.  For 
example,  should  the  ion  of  interest  become  surrounded  say 
by  more  electronegative  neighbors  the  valence  charge  would 
be  diminished  and  the  positive  potential  due  to  this  charge 
would  be  decreased,  effectively  lowering  the  entire  nuclear 
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FIG.  9  POTENTIAL  OF  CORE  ELECTION  AS  MODIFIED  BY  "SPHERICAL  SHELL"  VALENCE  ELECTRON 


33 


coulomb  potential  curve  seen  by  the  inner  electron  and 
increasing  its  binding  energy.  This  sensitivity  of  the  inner- 
core-electron  binding  energy  to  the  disposition  of  charge 
of  the  outer  valence  charge  is  termed  the  chemical  shift. 

The  foregoing  simplified  model  would  seem  to  indicate 
that  a  given  change  in  valence  charge  magnitude  would 
always  lead  to  the  same  chemical  shift.  Such  is  not  the 
case  in  reality;  the  valence  charge  is  not  removed  to 
infinity  as  assumed  above  but  Is  redistributed  amongst 
neighboring  ions  in  accordance  with  the  particular  nature  of 
the  ionic-valence  bonding.  There  is  also  an  interaction 
energy  due  to  some  wave  function  overlap  of  the  valence  and 
inner  core  electrons  which  changes  as  the  valence  charge 
redistributes.  Also  the  valence  electron  in  general  dees 
not  have  spherical  symmetry. 

There  is  therefore  a  chemical  shift  in  going  from  a 
free  atom  to  the  liquid  or  solid  state.  There  are  also 
chemical  shifts  when  ionic  environment  in  a  crystal  is 
changed  by  vacancies,  impurities,  or  order-disorder 
phase  transitions.  From  the  chemical  shift  one  may  try 
to  deduce  the  location  and  identity  of  surrounding  ions 
and  correlate  this  structural  information  with  device 
figures  of  merit  and  performance. 

c.  Features  of  APS.  APS  is  the  threshold  measurement  for 
electronic  excitation  of  an  inner  core  electron  of  a  surface 
atom  to  the  Fermi  level.  This  threshold  is  sharp  when  both 
the  inner  and  Fermi  levels  are  well-defined.  At  threshold 
both  the  incident  and  the  inner  core  electron  end  up  at 
the  Fermi  level  (leaving  a  localized  core  hole).  Since 
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this  transition  involves  the  disposition  of  two  electrons,  the 
excitation  probability  depends  on  the  square  of  final  state 
density,  i.e.  at  the  Fermi  level. 

Plasmon-aided  de-excitation  can  also  occur.  This  produces 
a  plasmon  satellite  structure  in  the  output  data,  with 
separations  determined  by  the  plasmon  energy.  Plasmon 
satellites  are  especially  strong  in  low  Z  materials  where  the 
core  size  is  such  as  to  couple  better  with  the 
electron  density  fluctuations. 

In  APS  the  experimental  procedure  probes  only  surface 
rather  than  bulk  atomic  phenomena  for  the  following  reason: 

The  moderate  energy  electrons  used  (200  to  2,000  eV)  do 
penetrate  many  monolayers.  But  the  incident  electron 
loses  some  energy  traversing  each  layer  and,  if  at 
threshold  energy  at  the  surface,  is  below  threshold  after 
penetrating  a  few  atomic  layers. 

d.  Experimental  Setup.  The  experimental  apparatus  is 
shown  schematically  in  Fig.  10.  Accelerated  electrons  from 
a  heated  tungsten  filament,  F,  strike  the  anode,  S,  which 
serves  as  the  sample  under  study.  As  the  anode  potential 
is  slowly  increased  there  is  a  certain  incident  electron 
energy  at  which  the  accelerated  filament  electrons  are 
energetically  capable  of  exciting  an  inner  core  electron  of  a 
surface  ion  to  the  Fermi  level;  at  this  threshold  both  the 
incident  electron  and  the  inner  core  electron  end  up  in  the 
first  unoccupied  states  above  the  Fermi  level. 

The  threshold  electron-excitation-level  for  a  core  level  of  a 
given  surface  ion  constitutes  the  primary  data  of  the  experiment. 


De-excitation,  i.e.  filling  of  the  core  vacancy  occurs 
both  by  radiative  and  Auger  processes.  Although  the 
latter  dominates  for  the  1  keV  energies  in  question  it  is 
nevertheless  the  gross  soft  x-ray  yield  which  is  measured 
in  APS;  complex  and  expensive  energy  analysers  for  decay  products 
are  thus  not  needed. 

The  soft  x-rays  pass  through  a  grid  in  front  of  a 
collection  can,  C,  and  generate  photoelectrons  upon  collision 
with  the  walls  of  the  can,  which  electrons  are,  in  turn, 
collected  by  a  wire,  W,  held  at  a  positive  potential. 

Auger  and  filament  electrons  are  kept  out  of  the  collector 
can  by  applying  a  negative  bias  voltage  to  the  can. 

The  experiment  is  conducted  as  follows:  As  the  sample 
voltage  is  increased  so  that  some  new  deeper  lying  core 
electron  energy  level  can  just  barely  be  excited  up  to  the 
Fermi  level,  the  sudden  "appearance"  of  a  small  increase  in 
x-ray  yield  is  a  measure  of  the  core  electron  energy  level. 

This  small  but  sudden  increase  is  more  easily  seen  amidst 
the  large  but  slowly  varying  soft-x-ray  background  yield 
(from  all  other  radiative  de-excitation  channels)  by 
resorting  to  electronic  differentiation  of  the  collector 
current.  This  is  done  in  the  usual  manner  by  superimposing 
a  "wiggle"  frequency  upon  the  sample  potential  and  using 
synchronous  detection  by  phase  lock  instrumentation. 

A  basic  limit  of  energy  resolution  is  the  energy  spread 
of  the  incident  electron  beam  due  to  the  voltage  drop  along 
the  filament  and  the  thermal  spread.  The  Auger  lifetime  of 
the  core  hole  also  provides  a  basic  linewidth. 
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An  APS  apparatus  can  be  rather  easily  modified  to  an 
Auger  Appearance  Potential  Spectroscopy  (AAPS)  apparatus 
by  adding  electron  collection  plates  nearby  the  sample. 

Although  energy  analysis  of  Auger  electrons  is  not  performed 
it  still  remains  a  "spectrometer"  in  that  inner-core 
energy  levels  can  be  determined.  Because  Auger  de-excitation 
dominates  in  the  soft  x-ray  energy  range  of  interest  these 
signals  would  be  more  easily  detected.  APS  and  AAPS  are 
complementary  techniques  in  the  sense  that  overlapping 
spectral  lines  from  two  different  surface  elements  in  say 
AAPS  would,  in  general,  be  well  separated  in  APS  and  vice  versa. 

An  ultrahigh  vacuum  system  has  been  built  and  a 
pressure  of  2.6  x  10-10  Torr  reached  after  a  200°C  bakeout. 

An  APS  spectrometer^^  has  been  built  on  a  Varian  vacuum 
flange  and  is  shown  in  Fig.  11.  The  electronic  voltage 
sweep  is  currently  being  designed  and  constructed.  The 
schedule  is  to  try  a  transition  metal  first  to  check  out 
spectrometer  calibration  and  sensitivity,  and  then  to 
proceed  to  the  lead  salts.  The  chemical  shift  for  say 
PbS  will  be  established  by  measuring  an  N  shell 
core  level  both  before  and  after  exposure  to  oxygen. 
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6*  Carrier  Scattering  In  Semiconductors  by  Screened  Surface 


Charges  (Theory:  R.  F.  Greene  and  John  Malamas#) 

Early  work  19  by  Greene  and  O’Donnell  In  19 66  had  established 
that  surface  charges,  e.g.  In  chemisorption  bonds,  provide 
the  dominant  surface  scattering  mechanism  at  fairly  well- 
prepared  semiconductor  surfaces.  Several  Important  problems 
remained  unsolved,  however,  until  the  present  work:1^ 

a.  The  actual  form  of  the  screened  potential  of  a 
partly  random  array  of  surface  charges  was  unknown  In  1966, 
so  that  an  arbitrary  Yukawa-type  of  scattering  potential  had 
to  be  used  then.  The  actual  three-dimensional  form  of  the 
space  charge  potential  of  such  an  array  was  determined 
finally  in  rn  earlier  part  of  this  program20  and  it  Is  this 
realistic  potential  that  is  used  in  the  present  wcrk.1^ 

b.  Dielectric  image  effects  were  omitted  in  the  Greene- 
0' Donnell  calculation.0^  This  force,  which  repels  carriers 
from  regions  of  lower  polarizability,  has  been  included  in 
the  present  calculation,  where  it  is  shown  to  decouple 
carriers  appreciably  from  surface  scattering  potentials. 

c.  Surface  charge  densities  are  commonly  high  enough 
for  several  scatterero  to  lie  within  an  electron  wavelength, 
resulting  in  strong  interferences  between  scattered  waves. 

This  was  ignored  in  the  Greene-0 ’ Donnell  calculation,  but 

is  treated  exactly  now.  This  is  done  by  treating  the  entire 
screened  array  as  a  single  scatterer,  whose  detailed  arrange¬ 
ment  is  represented  in  the  scattering  rate  by  the  structure 


*Army  Night  Vision  Laboratory  (not  supported  under  this  contract). 
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factor  S  of  the  array.  Numerical  values  of  the  Fuchs 
reflectivity  factor  are  presented,  for  various  1.  iterial 
parameters,  and  as  a  function  of  the  angular  coordinates 
of  electrons  leaving  the  surface.  Complete  details  of  this 
calculation  are  given  in  the  enclosure  (1),  which  is  a 
preprint  of  reference  16. 

7.  Surface  Transport  (Theory:  H.  F.  Greene) 

Surface  transport  calculations,  requiring  quantitative 
values  of  the  Fuch  s  carrier  reflectivity  and  its  angular 
dependence,  provide  values  of  surface  mobility,  surface 
Hall  coefficient,  and  surface  magnetoresistance,  with  which 
surface  scattering  models  can  be  tested  against  experiment. 
Improvements  are  needed  in  the  transport  theory,  principally 
because  the  standard  surface  transport  theory  does  not 
handle  the  energy  dependenc  of  bulk  scattering  time,  which 
is  crucial  for  the  surface  galvanomagnetic  properties.  We 
are  making  these  improvements  in  the  theory  and  carrying 
out  the  evaluation  of  the  surface  galvanomagnetic  coefficients 
numerically. 

8 .  Chemisorption  Thermodynamics  -  Surface  Scatt erlng  Connection 

(Theory:  R.  F.  Greene) 

We  have  made  an  important  breakthrough  in  the  theory  of 
semiconductor  surfaces,  with  the  discovery  of  a  direct 
relation  between  the  adsorption  isotherm  of  charged  chemisorbed 
atoms  and  the  surface  scattering  produced  in  the  semiconductor. 
(A  preliminary  form  of  this  result  has  been  published  in 
Thin  Solid  Films17  and  a  final  form  in  Surface  Science18). 

This  discovery  makes  it  possible  for  the  first  time  to  use 
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chemisorption  data  to  clarify  transport  phenomena  and  vice 

versa,  ana  should  greatly  increase  the  ability  of  scientists 
to  test  models  of  surface  structure. 

This  breakthrough  was  an  unexpected  result  of  oUr 
investigation  of  interference  effects  In  surface  scattering 
of  carriers  In  semiconductors.  Using  the  Born  approximation, 
we  found  that  these  interference  effects  could  be  expressed 
in  terms  of  the  structure  factor  S  which  represents  the 
arrangement  of  the  surface  scattorers.  We  used  the  fact 
that  S  Is  the  Fourier  transform  of  the  pair  correlation 
function  which  represents  the  effects  of  interaction 
between  adsorbed  atoms  in  the  statistical  thermodynamics  of 
chemisorption.  The  non-ideal  behavior  of  the  adsorbed  gas, 
as  measured  by  its  deviation  from  the  Langmuir  isotherm,  is 
directly  related  to  the  interference  effects  in  the  surface 
scattering  produced  by  these  adsorbed  atoms.  Both  the 
deviations  from  the  Lanmuir  isotherm  and  the  scattering 
interferences  are  obtainable  from  experiment,  and  a  test  of 
the  theory  will  be  attempted  experimentally  with  various  gases 
on  IV- VI  semiconductor  surfaces.  The  new  theoretical 
relation  can  be  written  in  terms  of  the  Fuchs  reflectivity 
P  (which  determines  the  surface  mobility)  as  follows: 

^  ^random  (kT/n(l-n) ) ( Sn/Sg)^ 

whe^  e  Prancjorn  is  the  calculated  value  of  p  for  a  random 
array,  n  is  the  mean  occupation  number  of  a  surface  site 
(coverage)  and  g  is  adion  chemical  potential,  obtainable  from 
the  overpressure.  Complete  details  are  given  in  enclosure  (2) 
which  is  a  preprint  of  reference  18. 
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?*  Theory  of  Chemisorption  Isotherms  on  Semiconductors 
(Theory:  R.  P.  Greene  and  C.  Richmond*) 


Work  has  begun  on  the  statistical  thermodynamic  calculation 
of  chemisorption  isotherms  for  various  gases  on  semiconductors, 
emphasizing  the  interactions  between  adions,  and  the  formation 
of  space-charge  layers.  Instead  of  the  familiar  adsorption 

isotherm  of  Langmuir  for  the  case  of  no  interaction  between  adions, 
which  can  be  written  in  the  form 


n  =  (1  +  exp ( ( eQ  -  g)/kT)-1 

(random  arrangement),  we  show  that  chemisorption  with  charge  transfer 
produces  a  non-Langmuir  isotherm  of  the  form 

H  =  (1  +  exp((eo  -  g  +  !e<J>p5(H))/kT)-] 

where  <pps  is  the  band-bending  potential,  with  which  we  express 
the  adion  interaction.  Numerical  evaluation  shows  that  this 
extra  electrostatic  interaction  depresses  the  coverage  n. 

Extensions  of  the  theory  are  being  made  to  take  into  account 
partial  ionization  of  the  chemisorbed  atoms,  and  to  treat 
the  interaction  by  the  Bragg-Williams  and  Bethe-Peierls 
methods. 

10*  D?-electric  image  Force  Theory  (Theory:  R.  F.  Greene) 

Cur  scattering  theory  has  shown  that  the  dielectric 
image  force  plays  an  important  role  in  surface  scattering. 

The  dielectric  image  force  is  also  known  to  be  important  in 

•Not  supported  under  this  contract. 
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heterojunction  tunnelling  and  in  surface  recombination  effects. 
Nevertheless,  the  theory  remains  in  an  unsatisfactory  state. 

The  present  classical  theory  diverges  as  the  electron  approaches 
the  surface,  and  the  sign  of  the  force  itself  is  unknown 
for  holes. 

We  have  begun  a  quantum  treatment  of  the  dielectric  image 
force,  introducing  dispersive  lattice  optical  modes 
explicitly  and  so  removing  the  classical  divergence.  Many 
body  techniques  are  being  employed  to  treat  the  nearly 
filled  band  case. 
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Scattering  of  Carriers  in  Semiconductors  by 
Screened  Surface  Charges 
R.  F.  Greene* 

Naval  Ordnance  Laboratory,  White  Oak,  Maryland  20910 

and 

John  Malamast 

Night  Vision  Laboratory,  Fort  Belvoir,  Virginia  22060 

ABSTRACT 

A  new  theory  is  given  of  surface  scattering  in 
semiconductors  by  surface  point  charges,  correspond¬ 
ing  to  trapping  in  chemisorption  bonds  or  defects. 
Interference  effects,  between  waves  scattered  from 
different  scatterers,  are  evaluated  in  terms  of 
the  statistical  structure  factor  of  the  surface 
point  charge  array.  Dielectric  image  effects  are 
included.  The  scattering  potential  is  taken  from 
the  recent  statistical  treatment  of  the  three 
dimensional  non-planar  semiconductor  surface  space 
charge,  given  by  Greene,  Bixler,  and  Lee.  The 
Need  for  multiple  scattering  treatment  of  evanescent 
states  is  discussed.  Comparison  with  previous  theory 
is  made. 

Enclosure  (1) 
(reference  15) 


INTRODUCTION 


In  .1  966  Greene  and  O’Donnell  proposed^"  that  the  apparent 
ffusenes s  of  semiconductor  surfaces  in  transport  measurements 
was  due  to  scattering  by  surface  charges,  e.g.  localized  at 
chemisorbed  atoms,  defects,  etc.  They  gave  a  simple  calculation 
of  this  scattering,  finding  a  rather  strong  angular  dependence. 

This  angular  dependence  seemed  to  explain  (Greene‘S)  the  lack  of  a 
surface  mobility  cusp  (Greene,  Frankl,  and  Zemel“’)  in  the  trans¬ 
port  measurements  of  Davis^  on  InSb  surfaces.  Further  confirmation 
of  the  GO  (Greene -0 ' Donnell)  theory  was  provided  by  the  measure¬ 
ments  of  Kamins  and  MacDonald^  on  the  Si-Si02  interface.  More 
recent  work  by  Preuss,6  however,  does  show  evidence  of  the 
surface  mobility  cusp,  suggesting  that  a  closer  look  at  surface 
scattering  is  needed.  Extensive  studies  of  scattering  at 
Si-Si02  interfaces  have  recently  been  reported  (e.g.  Sah7)  ,  also. 
Motivation  for  this  from  another  quarter  may  be  cited:  Recently 
it  has  been  shown  (Greene^)  that  there  exists  a  general  relation 
between  surface  scattering  rates  and  the  thermodynamics  of 
chemisorption,  where  the  chemisorption  bonds  provide  the  trapped 
charge  which  causes  the  scattering. 

The  GO  theory  is  oversimplified  in  several  respects:  a)  It 
neglects  the  interference  terms  in  the  scattering  from  different 
surface  charges;  b)  It  uses  an  arbitrary  form  for  the  screened 
potential  of  the  surface  charges;  and  c)  It  neglects  the  dielectric 
image  force.  The  present  paper  provides  a  new  theory  of  surface 
scattering  by  discrete  surface  charges  which  is  free  of  these 
limitations,  he  regard  the  entire  screened  array  of  surface 
charges  as  a  single  scatterer.  Interference  terms  can  then  be 
treated  exactly,  in  terms  of  the  statistical  structure  factor  of 
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the  surface  charge  array.  At  the  same  time  we  replace  the 

arbitrary  GO  scattering  potential  by  the  potential  calculated 

0 

by  Greene,  Bixler,  and  Lee,  (denoted  GBL  in  this  paper) 
of  the  inhomogeneous  surface  space  charge  region.  We  also 
take  into  account  the  GBL  dielectric  image  force  which,  in  its 
screened  form,  strongly  decouples  the  scattering  potential  from 
electrons  with  small  normal  momenta.  The  qualitative  features 
of  the  GO  theory  reappear,  but  with  somewhat  different  angular 
dependence . 

The  calculation  of  scattering  by  a  partially  ordered  array 
of  scatterers  becomes  tractable  only  with  use  of  the  statistical 
properties  of  the  scattering  potential.  This  is  also  supplied 
by  the  GBL  theory  which  treats  the  semiconductor  surface  space 
charge  potential  as  a  partially  correlated  stochastic  function 
of  position,  determined  by  the  surface  distribution  of  point 
charges  . 

The  density  of  surface  scatterers  is  often  high  enough  for 
scattering  interference  effects  to  be  very  strong.  The  concept 
of  the  differential  scattering  cross  section  of  individual 
scatterers  then  becomes  inexact  and  should  be  replaced  by  the 
dimensionless  differential  scattering  probability  of  the  surface 
as  a  whole. 

We  show  that  in  the  Born  approximation  all  interference 
effects  between  the  scattered  and  incident  waves  disappear,  for 
a  stochastic  scattering  potential,  but  the  interferences 
between  different  scatterers  remain.  The  need  for  a  multiple 
scattering  treatment  is  pointed  out,  particularly  for  the 
evanescent  or  surface  currents. 
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1.  Surface  Scattering  Theory  in  the  Born  Approximation 
We  represent  an  ideal  crystal  surface  by  an  infinite  potential 
step  at  z  =  0.  The  unscattered  electrons  are  then  described  by 


H%°  =  Et|)°  ;  z  >  0  ;  H°  =  -  fiV/2, 

\p°  =  (sin  k  z)e1^*^  ;  z  >  0 


(1.1a) 
(1. lb) 


<T  =  0 


;  z  <  0 


(1.1c) 


where  E  fi  kQ  /*.m,  and  kQ2  =  K2  +  k^2.  (We  write  2 -dimens ional 
vectors  as  capitals,  e.g.  K  =  (K  ,  K  ,  0).)  Next  we  introduce 

a  y 

a  scattering  potential  -e<J>(R,z)  which  produces  a  scattered  wave  \p: 


(H°  -  E  -  e  d) )  (ip°  +  ip)  =  0 ;  z  > 


ib  =  0 


;  z  <  0 


In  the  Born  approximation 


(H°  -  E)tP  =  e<^° 


;  z  >  0 


(1.2a) 

(1.2b) 


Cl. 3) 


It  is  convenient  to  use  the  2 -dimens ional  transforms  of  <J>  and  ip 

4>C<?,z)  =  / Cd2R/27T)  <J>(R,z)e'iQ'* 

4»CQ,  z)  =  y*Cd2R/2ir)  i|;(R,z)e"iQ*^ 

and  to  transform  Eq.  (1.3)  to  (Q,z)  space 


(1.4a) 


Cl. 4b) 


O  O  A  _ 

C3  /az  +  r  )^(Q,z)  =  ( -  2me/n  )  (sin  k,  z)  <p  (Q-K,z) ;  z  >  0  (1.5a) 


where 


1  “  C^Q2 -Q2) 17 2  for  kQ2  >  Q2,  and  T  =  i  (Q2-k  2)1/2  for  k  2  <  Q2 .  (1.5b) 


The  solution  of  Eq.  (1.5a),  subject,  to  Eq .  (1.2b)  can  be  written 


W~U“J  (  )  } 


^(Q,  O  =  (-2me/ii  )  I  dz'  G(z;z')  ^C^'K.z’)  sin  z  *  ;  z  >_  0.  (1.6a) 


Here  G(z,z')  is  the  one -dimensional  Green  function 


G(z,z')  =  (l/2iT)  (exp  if  jz-z'[  -  exp  ir|z+z'|  ) 


(1.6b) 


which  has  the  properties 


G(z,z')  =  0  for  z  =  0 


(1.7a) 


and  (32/3z2  +  r2)  G(z.z’)  =  <5(z  -  z')  +  6 (z  +  z'), 


(1.7b) 


which  ensure  that  both  Eq.  (1.2b)  and  (1.5a)  are  satisified. 

The  choice  of  (+)  exponents  in  G(z,z')  means  that  we  are 

2  2 

requiring  that  the  scattered  propagating  waves,  Q  <  k^  , 


all  have  the  outgoing  form  e 


+  irz 


for  z->°°,  as  seems  physically 


reasonable.  Moreover,  were  we  to  admit  incoming  scattered  waves, 

2  2 

these  would  produce  divergent  exponential  solutions  for  Q  >  kQ  . 

The  GBL  potential  produced  by  surface  point  charges  consists  of 
a  planar  or  band-bending  term  which  depends  on  z  only,  and  an 
inhomogeneity  term,  <J>(R,z),  which  produces  scattering  and  has  the 
form 


<KQ,z)  =  <KQ,0)  e , 


where 


<  =  (q 


2  +  q2)1/2, 


(1.8a) 


(1.8b) 


and  where  q  is  the  effective  inverse  screening  length  near 
the  surface: 


qs2  =  (47re2nB/e  +  kT)  exp(e<J>  /kT)  =  qQ2  (exp  e4>  /kT) 


(1.8c) 


x.  -  —  *  u-S  (  |  i 


is  the  value 


Here  nn  is  the  bulk  electron  concentration  and  d> 

■D  ps 

of  the  band-bending  potential  at  the  surface.  It  might  be 
objected  that  the  use  of  the  Fermi -Thomas  approximation  in  GBL  theory 
is  inconsistent  with  our  scattering  model,  viz.  Eqs.  (1.1c)  and 
(1.2b),  whereby  carrier  wave  functions  are  terminated  at  the  surface. 
It  has  been  shown  (Greene1  )/the  screening  approximations  are  neg- 
igable  for  surface  fields  below  about  5  x  10^  V/cm  in  non¬ 
degenerate  semiconductors.  If  we  use  the  simple  z-dependence 
of  Eq.  (1.8a),  then  we  cm  write  the  scattered  solution  Eq.  (1.6a) 
in  the  explicit  form 

K3,z)  =  (me/2fi2)  <KQ-K,0)  (Fp  eirz+  Fg(z)),  (1.9a) 

where  F_  =  F*  =  8<k  /(£k2  +T2  +  k  2)2  -  4k  2r2)  ,  (1.9b) 

pp  z  z  z 

r  i(k _-<)z  -  2  2 

and  F  =  F*  =  2i  e  2  /(<  -k/  +  TZ  -  2i<k  )  +  c.c.  . 

SSI  Z  Z  j 

(1.9c) 

Here  k  =  k(Q-K)  =  ((Q-K)2  +  qs2)1/2  (1.9d) 

and  T  is  given,  by  Eq.  (1.5b).  The  Fp  term  is  a  propagating 
wave  for  Q2  <  k  2,  but  is  evanescent  for  Q2  >  k  2,  The  F 
term  is  evanescent  for  all  values  of  Q.  Electrons  scattered 
into  evanescent  states  clearly  interact  strongly  with  the 
scattering  potential  and  will  be  multipll -scattered.  A 
treatment  of  these  states  will  be  given  elsewhere. 
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2.  Scattered  Flux  Density  from  a  Partially  Random  Scatterer 
We  now  calculate  the  scattered  flux  density  J,  and  then 
afterwards  introduce  the  statistical  or  random  properties  of 
the  scattering  potential.  The  probability  flux  density  for 
a  state  rp  =  \p°  +  \pt  viz 

J  =  (ifi/2m)  (TVT*  -  c.c.),  (2.1) 

has  three  distinct  terms 

j  =  J C^P° ;  t/»°)  +  +  J(i|>°;*jO  (2.2) 

corresponding  to  the  unscattered  wave,  the  scattered  wave, 
and  the  interference  between  the  scattered  and  unscattered 
waves,  respectively.  These  take  on  complicated  forms  when 
the  scattering  potential  has  the  compli  cated  form  of  the 
inhomogeneous  space  charge  potential  $(R,z) . 

Simplifications  occur,  however,  when  we  make  use  <H  the 

statistical  properties  of  the  scattering  potential.  For  our 

-+ 

purposes  <j>(R,z)  may  be  treated  as  a  random  function  of  R  with 
two  basic  ensemble  average  properties: 

(<KR,z)}  =  0  ,  (2.3a) 

and  (<j>(^,z)  (R+  AR ,  z )  =  function  of_^AI$  and  z 

but  not  of  R.  (2 . 3b) 

Eq.  (2.3a)  indicates  that  we  have  subtracted  out  the  mean  value 
of  the  space  charge  potential,  as  mentioned  above  in  connection 
with  Eq.  (1.2a).  Eq.  (2.3b)  means  that  the  various  probability 
density  functions  for  (R)  are  uniform  in  R.  Because  of  thi  - 
statistical  uniformity,  we  may  equate  ensemble  averages  with 
spatial  averages,  e.g. 


/ 


{ <f>  (R) }  =  4>  =  lim  L 

L-^oo 

For  use  in  our  scattering  theory  it  is  convenient  to  re-express 
conditions  (2.3a)  and  2.3b)  in  Q-space.  This  is  possible 
because  each  member  of  the  statistical  ensemble  of  functions 
4>(R,z)  can  be  represented  just  as  well  by  its  Fourier  transform 
<KQ» z) .  Equations  (2.3a)  and  (2.21;)  then  appear  as 

{ <J>  (Q ,  z) }  =  0  (2.5a) 


/ 


d^R  cf>  (R) 


(2.4) 


and 


(<f>  (Q » z)  (Q  '  ,  z)  *}  =  4tt  6  CQ  -  Q’)  g^CQ.z), 


(2.5b) 


where  g^(Q)  is  the  spectral  density  of  tj>(R,z) 


-  O 


g,(Q,z)  a  lim  L  j  <{>T  (Q  ,  z) 


and  H  0  for  R2  >  L2  (2.5c) 


2.5a)  follows  directly  from  ensemble-averaging  Eq.  (1.4a)  and  using 
Eq.  (2.3a).  Eq.  (2.5b)  can  be  obtained  by  writing 

(<KQ)  4>  (Q  '  )  * )  =  j (d2R/ 2tt)  J  (d2R'/27r)  {^(R)0(R')f  exp(iQ-R  -  iQ'-R') 
and  then  imposing  the  statistical  uniformity  condition  (2.3b),  giving 

U(Q)<KQ')*>  =  <5  f  Q  -  Q')  f  (d2AR/2Tr)  p<jj(AR)  exp  -  iQ-hl. 

where  p.(AR)  =  lim  L  2  /  d^R  cf>L(R,z)  4>.  (R+AR,z) 

Eq.  (2.5b)  then  follows  by  means  of  the  Wiener-Khinchine  relation, 

(Ming  Chen  Wang  and  Uhlenbeck^) . 

We  can  now  use  Eqs.  (2.5a)  and  (2.5b),  the  statistical 

properties  of  the  space  charge  potential,  to  simplify  the 

probability  flux  density  J,  when  the  latter  is  ensemble  -  averaged . 

First  we  note  that  in  the  Born  approximation  '<p  is  a  linear 

functional  of  so  that  . 

. -1  (  /  /. 


8 


(2.6) 


I 


=  0 


Particle  flux  conservation  in  surface  scattering  requires  that 

V^}  vanish  exactly,  since  carriers  are  assumed  neither 
to  go  through  the  surface  nor  to  be  created  there.  Of  course 


W  =  0  because  \pQ  is  a  standing  wr 


Furthermore,  the 


exact  i|>  contains  only  outgoing  (k  >  0)  Waves  for  z  -*■  00  so  that 
(J  (i|>,\|0)  >  0.  Therefore  the  exact  interference  terms  { Jz (ip  ,ipQ) }  +  c.c 
must  be  negative.  Eq.  (2.6)  then  means  that  our  approximate  ip  fails 
to  conserve  particle  flux:  this  will  be  discussed  elsewhere.  Next, 
we  point  out  the  great  simplification  produced  by  Eq.  (2.5b)  in  the 
scattered  flux  density  J(^,ip)  whereby 


J(iM)  =  (ih/2m)  (me/2h)  2  j  (d2Q/27r)  J  (d2Q  '  /  2tt)  d>  (Q- -K  ,  0)  4.  (Q '  -1 , 0)  * 


[((FTieirz  +  Fs)ei^^)V((F!e'ir*z  +  Fi)e'iQ,’R)-c.c.] 
P  ^  P  ^ 


(2.7) 


is  reduced  to  an  expression  which,  as  shall  show  below,  has 
a  simple  semi-classical  interpretation  and  is  capable  of 
complete  evaluation: 

UOM)}  =  (me2/8h2)  J  d2Q  g^ (Q-K,0) [  (2Q+lz  (T+r*) )  F2  exp  (ifz -ir*z) 


+  2Q  F2+2QFpFs(eirz+c.c)  -  lz  FpFse'KZ] 


where 

\  V  n 

*  1  V  7  7  ' 

-F  =  [(<-ik _)(e1  -c.c)  +  i(e  -c.c)]M  e  “  -c.c 

s  z 


(2.8a) 


(2.8’d) 


M  =  2i ( (k- ikz) 2  +  T2)'1,  and  tz  =  (0,0,1) 


(2.8c) 
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The  Fp  term  is  propagating  for  Q2  <  k2 ,  evanescent  for  Q2  >  k2. 
The  Fg  and  Fg  terms  are  evanescent  for  all  values  of  Q.  These 
evanescent  fluxes,  corresponding  to  real  surface  currents, 
are  not  correctly  given  in  the  Born  approximation  because  the 
re-scattering  of  electrons  in  such  states  cannot  be  neglected 
in  any  approximation.  In  this  paper  we  will  therefore  treat 
only  the  propagating  flux  density 


{J}p  =  (me2/4fi3) 


( 


Q2<k2 
x  0 


d2Q  f2  gtCQ-K,0HQx,  Qy,  r) 


(2.9) 


This  takes  on  a  clear  semi-classical  form  in  terms  of  the 
angular  coordinates  of  the  scattered  waves 


Qx  =  kQ  sin  8  cos  p  ,  Qy  -  kQ  sin  6  sin  n  . 


Writing  d2Q 
one  gets 


Q  dQ  dp 


ydy  dp, 


(y  =  cos  8) 


1 


f(K,-kz 


;y ,n) v 


where 

v  =  (Ti/m)CQx,  Qy,  T) 
is  the  scattered  velocity  vector,  and 


(2.10a) 

(2.10b) 


(2.11a) 


(2.11c) 


f(K,-kz;y,p)  =  (me/2h2)2yk^  F2  g^(Q-K,0)  (2.11c) 

is  the  scattered  distribution  function,  the  number  of  classical 
electrons  per  unit  solid  angle  at  (y,p). 

Although  it  has  a  clear  semi -classical  interpretation,  Eq.  (2.11a) 
is  somewhat  more  complicated  than  the  corresponding  expression  in 
the  theory  of  scattering  from  a  finite  range  potential  (Schiff12) . 


-- -  (  I  ) 
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There,  only  one  propagation  vector  is  seen  at  a  distant  detector  because 
the  scatterer  subtends  only  an  infinitesimal  solid  angle  at  the 
detector.  The  range  (2tt)  of  scattering  angles  appeals  in  Eq. 

(2.11a)  because  the  scattering  surface  subtends  the  solid  angle 
(2u)  at  a  detector,  no  matter  how  far  from  the  surface  it  may  be. 


3.  Differential  Scattering  Probability 
The  concept  of  differential  scattering  cross  section,  appropriate 
for  finite  range  scatterers,  is  generally  not  useful  for  surface 
scattering  except  when  the  surface  density  of  scatterers  is  low 
enough  to  ignore  interference  effects  and  screening  interactions. 

The  generally  appropriate  concept  is  that  of  the  dimensionless 
differential  surface  scattering  probability  w(3!, -k  ;y,p)  which  we 
define  as 

w(K,-kz;u,n)  =  N(y,n)/No(TE,-kz)  (3.1) 


where  NQ(K,-kz)  is  the  number/sec.  of  electrons  incident  on  unit 

area  of  the  surface  with  momentum  (Kv,  K„,  -k.)  ,  and  N(y,n)dydn 
_____  x  y  z 

is  the  number/sec.  of  particles  scattered  into  solid  angle  dydn 
from  unit  area  of  the  s  urf ace . 

Equation  (2.11a)  enables  us  to  find  N(y,n)  directly.  Evidently 

N(y,n)  =  y  v  f (y,n)  •  (3.2) 

On  the  other  hand,  the  unperturbed  wave  function  ip°  of  Eq.  (1.1b) 
describes  an  incident  flux  density  of  magnitude  v/4,  so  that 

V*,-kz)  =  Ukv/4,  where  y£  =  -  kz/kQ  <  0  .  (3.3) 

We  thus  obtain 

w(K,-kz;u,n)  =  (-y/yjj)  f  (y,n)  =  (-y  2/h^)  (mekQ/h2)  g^(Q-K,0)  ,  (3.4) 

where  is  given  by  Eq.  (1.9b)  and  g^(Q-K,0)  is  the  spectral  density 
of  the  scattering  potential  for  the  surface  momentum  change  (Q  -  K) . 
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This  spectral  density  g^CQ.z)  is  expressed  in  the  GBL  theory 
in  terms  of  the  spectral  density  g^(0)  of  the  surface  charge  density 
or  rather  of  <$ft  =  ft(R)  -  ft  (by  means  of  Poisson's  equation, 
the  static  electron  response  function,  and  the  Maxwell  boundary 
conditions) 


g^c^.z)  =  #(Q)2  e'2,cz  gM(Q)  (3 . 5a) 

where 

ft($)  =  (4-n/ e  +  )  2  (<  +  yQ)’2,  Cy  =  e  _/e  +  <  1)  (3.5b) 

describes  the  response  of  the  semiconductor  surface  space  charge 
region  to  an  arbitrary  surface  charge  distribution.  Suppose, 
in  particular,  that  ft(R)  consists  of  an  array  of  point  charges, 
e.g.  associated  with  chemisorption  bonds,  occupying  some  of  the 
sites  A  of  a  regular  surface  lattice  (mesh).  In  that  case  gr^ 
is  essentially  the  statistical  structure  facto’  S(Q)  of  that 
array: 

g^lQ)  =  (a  ^2/4Tf2)  S(Q)  (3.6a) 

S(Q)  =  E  (p(AA)/n  -  l)e1^'A^  +  E  6  ((J-S. )  (4Tr2/a)  (3.6b) 

At  Qi^O 

Here  n(A)  =  0,1  is  the  occupation  number  of  site  A,  n  is  the  mean  of 
n  (A)  , 


n  =  (n)  =  aft  , 


(3.6c) 


where  a  is  the  unit  mesh  area. 


Also,  p(AA)  is  the  pair  correlation 


function 


P (AA)  =  (n (A)n (A+AA) }  /n  ;  p(0)  S  1  . 


(3.6d,e) 


and  are  the  reciprocal  mesh  vectors 

Q^* A  =  2 it  integer,  or  zero.  (3.6f) 

These  reciprocal  mesh  vectors  determine  the  features  of  the 
evanescent  states  but  do  not  directly  appear  in  the  propagating 
flux  {J)p,  and  so  can  be  omitted  in  this  paper.  We  thus  obtain 
the  explicit  form  of  the  differential  surface  scattering  probability 

2 

W(p,n)  =  wG0(y,n)  ( - - -A  n  S(Q-K)  (3.7a) 

U  +  Y  |Q  -  R | / 


(3.7b) 


is  the  corresponding  result  of  the  GO  calculation.  The  factors 
multiplying  wqq  express  the  effect  of  taking  into  account  properly 
the  screening  potential  and  the  interference  of  waves  scattered 
from  different  scatterers. 

The  simplest  interference  situation  to  analyze  is,  of  course, 
that  of  no  correlation  at  all  in  the  occupation  numbers:  If 


then 


p(AA)  =  h  for  all  &R  ?  0 
S°(Q)  =  S° (0)  =  (1  -  n)/n  . 


(3.8a) 

(3.8b) 


14 


ihen  if  n  «  1,  w  is  linear  in  n  and  an  effective  surface  scattering 
cros_  section  might  be  introduced.  Next  consider  the  effect  of  a 


short  range  correlation:  If 

p(AA)  =  ii  for  all  |A^|kQ  > 

1/10 

(3.9a) 

then 

S (Q)  =  S(0)  =  X)  CpCAA)/n  - 

1)  • 

(3.9b) 

AA 


In  both  of  these  cases  the  scattering  is  determined  by  S(0)  which, 
as  has  been  shown  elsewhere  (Greene  )  is  a  thermodynamic  property 
of  the  chemisorbed  charged  adatoms,  if  these  have  reached  equilibrium. 
If  the  distribution  has  not  reached  equilibrium,  so  that  S(0) 
cannot  be  obtained  from  statistical  thermodynamic  arguments, 
some  guidance  in  the  choice  of  p(A^)  is  available  from  the  sum  rule 

fd2Q  S(3)  =  (4u2/ii  a)  (1  -  n)  (3.10) 

BZ 

in  which  the  integral  goes  over  one  Brillouin  zone  of  the 
reciprocal  mesh,  and  which  arises  from  Eq.  (3.6e). 

Finally,  we  consider  the  space  charge  form  factor, 

B  =  4k2  (k  +T]Q  -  K|;T2,  (3.11) 

appearing  in  Eq.  (3.7a),  which  corrects  for  the  arbitrary  potential 
used  in  the  GO  treatment,  viz  (e/e  +  r)  exp(-qQr),  for  z  >_  0. 

For  small  angle  scattering  events  one  can  see  that 

B  =  4,  for  |Q  -  K |  <<  qg  (3.12) 
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U  ) 


Thus  the  Fourier  coefficients  of  the  Yukawa  potential  mentioned 

above  are  exactly  half  the  correct  ones,  for  Q  ■+  0.  At  the  other 

2  2 

extreme,  one  can  consider  large  angle  scattering  for  qs/kQ  <<:  1> 
in  which  case 

B  =  4(1  +  y)  2  for  |Q  -  K|  »  qg .  (3.13) 

This  simply  menns  that  large  angle  scattering  events  are  determined 
by  the  unscreened  Coulomb  singularity  itself,  which  has  the  form 
(e/c+r)  (2/(1  +  Y)) . 


-  (  /  ) 
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4.  Dielectric  Image  Effects 

The  polarization  energy  of  an  electron  changes  as  the  electron 
is  brought  nearer  the  interface  between  two  dielectrics.  This 

1  7 

can  be  expressed  (Abraham  and  Becker  )  in  terms  of  a  dielectric 
image  force,  and  has  already  entered  the  scattering  potential  <f> 
directly,  as  is  shown  by  the  factor  y  =  e_/e+  in  the  spectral  density 
g^  in  Eq.  (3. 5b)  . 

But  there  is  also  a  direct  dielectric  image  repulsion  of  each 
electron  from  the  surface  when  y  <  1.  In  GBL  it  was  shown  that  tor 
semiconductors  the  work  done  on  an  electron  in  bringing  it  and  its 
screening  cloud  closer  to  a  surface  can  be  expressed  in  terms  of 
a  screened  image  potential  <j>j  given  by 

-eifrj  ='  (b  e2/4e  +  z)  exp  -qgz  ,  (4.1) 

where  b  is  a  constant;  between  unity  and  (y  -  l)/(y  +  1). 

Being  independent  of  R,  this  is  not  a  scattering  potential. 
Nevertheless,  the  scattering  is  weakened  because  <J>j  decreases  the 
amplitude  of  the  unscattered  wave  functions  ij)0  near  the  surface. 

This  direct  dielectric  image  repulsion  can  be  introduced 
into  the  foregoing  scattering  theory  very  simply,  albeit  crudely, 
by  terminating  at  the  classical  turning  point  Zj  instead  of 
at  the  surface  z  =  0,  z^  being  given  by 

-  e<j>j(Zj)  =  ft2k2/2m  (Zj  >  0)  .  (4.2) 

Electrons  incident  at  grazing  angles  have  small  kz  and  hence  have 
large  values  of  Zj  and  are  strongly  de-coupled  from  the  scattering 
potential.  We  thus  replace  Eq.  (1.1)  by 
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. (  I  ) 


t|)°  =  (sin  kz(z  -  zI))elK*R,  z  >  Zj  (4.3a) 

ijj°  =  0,  z  <  z  j .  (4.3b) 

For  the  Born  approximation  treatment  of  Section  (1)  to  remain  valid, 
we  must  correspondingly  replace  Eq.  (1.2b)  by 

\p  =  0,  z  <_  Zj  (4.4) 

The  Born  approximation  Eq.  (1.5a)  is  then  replace  by 

(32/3z2  +  r2)ip(Q,z)  =  ( -2me/h2)  sin  kz(z-Zj)  0(Q-K,z)  ,  z  >_  Zj.  (4.3) 
The  scattered  wave  Eq.  (1.6a)  is  now  replaced  by 

v*  00  ^ 

\p  (Q,  z)  =  (-2me/h2)  /  dz  •  GD(z  ,  z  • )  <p  (Q -K ,  z  ' )  sin  k^z'-Zj)  (4.6) 

ZI 

where  now  the  displaced  Green  function 

Gjjtz.z’)  =  (l/2iT)  (exp  iT |  (z-Zj)  -  (z ' -Zj)  | 

-  exp  ir|(z-Zj)  +  (z-Zj)|)  (4.7a) 

has  the  property 

G^(z,z')  =  0  for  z  =  Zj  (4.7b) 

(which  ensures  that  Eq.  (4.4)  is  satisfied),  and  the  property 

(32/3z2  +  r2)GD(z,z')  =  6 (z-z ' )  +  6  (z  +  z  ’  -2Zj)  (4.7c) 

(which  ensures  that  Eq..  (4.5)  is  satisfied.)  Here  the  fact  that 
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. 4(1), 


4>  (Q-K,  z  ' )  =  <f>  (Q  -  K ,  0 )  exp(-Kz')  =  4>CQ-K,Az)  exp(-tczj), 


(4.8) 


where  Azf  =  z'  -  Zj,  allows  Eq.  (4.6)  to  be  written  as 

Th  (Q,  Z)  =  (exp-KZT)  (2me/1i2)  f  dAz’  G(Az  ,Az  ’ )  sin(kzAz  ’)4>  (^-K,Az  ' ) 

JQ  (4.9) 

where  now  G(Az,Az’)  has  exactly  the  same  functional  form  as  the 
Green  function  used  in  Eq.  (1.6a)  for  the  scattered  wave  without 
this  dielectric  image  effect.  Thus  the  entire  effect  of  shifting 
the  unscattered  wave  a  distance  zj(^z)  as  merely  to  shift  all 
the  scattered  waves  by  the  same  amount  and  to  change  their 
amplitudes  by  exp  (-<Zj). 

Because  this  factor  is  independent  of  3  and  z,  the  scattered 
current  density  is  unchanged  except  for  the  appearance  of  a 
multiplicative  factor  exp(-2icZj).  Thus  for  electrons  the 
differential  surface  scattering  probability  of  Eq.  (3.7a)  becomes 

w(y  ,n)  -  wGQ(y,n)  B  ^GXP  ~2<zl^n  S(^-K)  (4  10) 

where  B,  the  space  charge  form  factor,  is  given  by  Eq.  (3.11). 

Here  cos-1y  and  n  are  the  angles  of  the  emerging  electron  in  spherical 

coordinates . 

As  noted  above,  the  dielect"ic  image  factor  (exp-2iczj  (k^)) 
strongly  decreases  the  scattering  rate  for  electrons  at  grazing 
incidence,  because  Zj  increases  as  kz  decreases.  Inasmuch  as 
K  =  (q2  +  (Q  K)2)1/2,  it  is  clear  that  this  suppression  of 
scattering  occurs  mostly  for  large  angle  scattering.  This  is 
indeed  plausible,  since  the  spatial  shift  of  the  unscattered 
wave  primarily  cuts  down  the  overlap  of  ip°  with  the  more 
singular  part  of  the  scattering  potential. 
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5.  Fuchs  Reflectivity;  Conclusions 
We  use  the  general  boundary  condition  formalism  (Greene‘S) 
to  get  the  Fuchs  reflectivity  p  and  kinetic  specularity  WQ  from  the 
differential  scattering  probability  of  Eq.  (4.10):  For  an 
isotropic  surface 


WQ(y+)  =  l  -  JduJ  dAn(-y_/y+)ws(y_,n_ |y+,nj 

-1  0 

r°  r2v 

p(y+)  =  1  -  I  dy  /  dAn(-y./y+)wsCy.,n. |y+,n+) [l 
-1  Jo 


co 'An 


(5.1a) 


’  1  ■  y  2 \  1/2 


r_T 

!y+ 


(5.1b) 


where  cos  ^y  is  the  angle  between  the  scattered  k-vector  and  the 
inward  normal  to  the  surface,  and  An  =  n_  -  n+  is  the  angle 
between  surface  components  of  scattered  and  incident  k-vectors. 
These  integrations,  which  could  be  carried  out  analytically 
in  the  simpler  GO  theory,  unfortunately  had  to  be  carried  out 
numerically  for  wg  given  by  Eq.  (4.10).  We  carried  out  the 
integrations  for  the  completely  uncorrelated  ("pure  random") 
distribution  of  surface  charges,  for  which  the  structure  factor 
(leaving  out  the  Bragg  terms  outride  the  Brillouin  zone)  is 


S°(Q)  -  S°(0)  =  (1  -  n)/n  ,  (for  p(AA)  =  n  for  AA  /  0)  (5.2a) 

ii.  which  case  Eq.  (4.10)  reduces  to  the  simpler  form 

v;  =  wGQ  B(exp  -  2<Zj)  (1  -  n)  (5.2b) 

in  which  case  1  -  W  and  1  -  p  are  linear  in  n  for  n  <<  1. 
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ENCLOSURE  (  |  ) 


The  scatterers  are  charged  so  that  an  urxorrelated  distribution 
is  not  really  to  be  expected  except  when  h  is  close  to  zero. 

Correlation  effects  can  depress  S(Q)  and  wg  very  significantly, 
as  we  now  show  for  the  case  of  depressed  nearest  neighbor  occupation 
p(l)  <  n  of  a  simple  square  mesh.  From  Eq.  (3.6b) 

S(Q)  =  (1  -  n  -  2(n-p(l))  (cos  QxAq  +  c,;s  QyAQ))/n  >  0.  (5.3) 

In  the  absence  of  a  statistical  thermodynamic  calculation  of  p(l), 
we  might  guess  a  form  (which  keeps  S(Q)  positive  definite). 

n  -  p(l)  =  (1/4) (1 -n) (1 -e"a)  (5.4a) 

where 

a  =  2e2/Ao(e+  +  eJkT  (5.4b) 

is,  roughly,  the  nearest  neighbor  Coulomb  energy  in  units  of  kT. 

Since  in  wg  only  small  values  of  Qx,  Qy  'v  kQ  <<  |Aq|  ^  occur,  we  can 
set 

S(Q)  =  S(0)  =  (1  -  n)  e-ct/n  =  Sor0)e'a  (5.5) 

so  that  our  calculated  (uncorrelated)  values  of  1  -  W  and  1  -  v 

s  r 

have  to  be  multiplied  by  the  very  strong  factor  e~a. 

The  calculations  have  been  carried  out  for  spherical  energy 
surfaces  and  so  are  most  directly  applicable  to  semiconductors  like 
InSb  (Willardson15)  whose  electron  effective  mass,  0.013mo,  and  dielectric 
constant,  e  =  17,  were  used.  Numerical  results  for  p  and  W 

o 

versus  the  emergence  angle  cosine  y+  are  plotted  in  Figures  1 
through  4  for  three  different  forms  of  differential  scattering 
probability  wg .  Curves  labelled  GME  correspond  to  the  use  of 
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Eq.  (4.10).  Curves  labelled  GM  correspond  to  the  use  of  Eq.  (4.10) 
with  the  dielectric  image  factor  (exp  -2kZj)  omitted,  i.e.  to 
Eq.  (3.7a).  Curves  labelled  GO  correspond  to  the  GO  theory,'1'  i.e. 
to  the  use  of  of  Eq.  (3.7b).  In  Fig.  5  we  show  the  specific 
effect  of  the  dielectric  image  factor,  expressed  as 

{exp  -  2 k z r }  =  (1  -  pGME)/(l  -  pGM),  (5.6) 

plotted  as  a  function  of  the  emergent  angle  cosine.  The  dielectric 
image  factor  is  more  effective  for  slower  particles  and  nearer 
grazing  emergence  because,  in  both  cases,  k„  is  smaller. 

The  qualitative  features  of  the  simpler  GO  theory  are  retained: 
a),  p  +  1  for  grazing  emergence;  and  b) ,  the  Fuchs  diffusivity  (1  -  p) 
is  much  smaller  than  the  total  scattering  probability  (1  -  W  ). 

But  the  use  of  a  realistic  1 reatment  of  the  screening  of  the  surfac: 
scattering  potential  has  increased  its  Fourier  transform  and  hence 
(1  -  p)  by  about  a  factor  of  four.  This  is  partially  compensated 
by  the  dielectric  image  effect,  which  decreases  the  Fuchs  diffusivity 
(1  -  p) ,  particularly  at  grazing  angles.  For  such  angles  kz  is 
small,  Zj  is  large  and  (exp  -  2kZj)  becomes  much  less  than  unity. 

Compared  with  the  GO  theory,  as  Figures  1  thru  4  show,  the 
present  theory  gives  a  weakened  angular  dependence  of  p  at  grazing 
angles.  This  makes  cuspid  behavior  in  the  surface  mobility 
somewhat  more  understandable r 

One  unwelcome  ieature  of  the  GO  theory  which  is  retained  in 
the  present  theory  is  the  breakdown  for  larger  values  of  n,  where 
p  can  become  negative.  We  can  now  see,  however,  that  this 
breakdown  is  not  due  to  the  neglect  of  interference  effects 
but  has  to  be  ascribed  to  the  use  of  the  Born  approximation  with 


( 
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its  neglect  of  multiple  scattering.  This  is  particularly  serious 
for  evanescent  states,  a  treatment  of  which  will  be  given  elsewhere. 

A  Boltzmann  transport  calculation  of  surface  mobility  using 
the  calculated  values  of  the  Fuchs  p  now  becomes  possible  and  is  being 
carried  out.  That  will  then  permit  the  scattering  model  of  the  present 
calculation  to  be  tested  against  surface  mobility  measurements, 
e.g.  those  of  Preuss  on  InSb,  those  of  Sail  on  the  SitSiC^  interface, 
etc.  It  should  also  be  possible  to  test  the  present  surface  scatter¬ 
ing  theory,  or  rather  the  surface  charge  structure  factor  deduced 
from  a  fit  to  transport  data,  against  structure  factor  information 

from  chemisorption  isotherm  data,  using  the  new  surface  scattering- 

0 

chemisorption  relation  of  one  of  us  (Greene  ) .  Such  experiments, 
involving  chemisorption  and  electrical  transport  studies  in  thin 
single  crystal  films  of  IV-VI  semiconductors  with  exposure  to  oxygen, 
hydrogen,  etc.  are  being  carried  out  at  Naval  Ordnance  Laboratory 
(Lee^6) .  Further  scattering  calculations  are  also  being  pursued 
replacing  the  Born  approximation  by  a  one-particle  Green  function 
method  which  conserves  particle  flux  exactly,  and  which  includes 
ellipsoidal  energy  surface  effects. 


•3  U  ) 
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FIGURE  CAPTIONS 


Fig. 

la: 

p  versus  y.;  n/a  =  10^  cm  ^ ,  fi^k^/2m*  =  0.026  eV. 
m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  105  cm'1. 

Fig. 

lb: 

WQ  versus  y+;  n/a  =  10^  cm  ^ ,  fi^k^/2m*  =  0.026  eV, 
m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  105  cm'1. 

Fig. 

2: 

p  and  WQ  versus  y+;  n/a  =  10^  cm  ^ ,  -h^k^/2m*  =  0.026  eV 
m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  10^  cm  1. 

Fig. 

3: 

p  versus  y+;  n/a  =  10^  cm  fi^k^/2m*  =  0.26  eV, 

m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  105  cm'1. 

Fig. 

4 : 

p  versus  y  ;  ii/a  =  10^  cm  ^ ,  1i^k^/2m*  =  0.0026  eV , 

'  *T 

m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  105  cm'1. 

Fig. 

5: 

2  2 

{exp  -  2kZj}  versus  y+  for  different  values  of  1i  k  /2m* 
m*/m  =  0.013,  e+  =  17,  qQ  =  1.33  x  105  cm'1. 

i  ::-gcux£ 
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The  Connection  Between  the  Thermodynamics  of 
Chemisorption  on  Semiconductor  Surfaces 
and  Surface  .Scattering  of  Carriers 

Richard  F.  Greene* 

Naval  Ordnance  Laboratory 
White  Oak,  Maryland  20910 

ABSTRACT 

An  explicit  theoretical  relation  is  derived 
which  connecto  the  thermodynamics .of  chemisorption 
on  semiconductor  surfaces  with  the  scattering  rate 
of  conduction  and  valence  band  carriers  by  charger 
localized  in  the  chemisorption  bonds. 
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(reference  17) 


I.  INTRODUCTION 


In  surface  science  it  is  desirable  to  make  connections 
between  "exterior  phenomena",  such  as  chemisorption,  and 
"interior  phenomena",  such  as  transport  in  the  semiconductor 
space  charge  layer.  Of  course,  it  is  well  known  that 
chemisorption  bonds  are  often  charged,  producing  surface 
space  charge  layers  of  appreciable  conductance,  *  and  it  is 
now  also  accepted  that  chemisorption  bond  charges  are  strong 
scatterers  of  these  carriers.  In  this  paper  we  derive  an 
explicit  theoretical  relation^  which  ties  together  chemisorption 
thermodynamics  and  surface  scattering  rates. 

We  have  recently  calculated"'’  the  differential  surface 
scattering  rate  w,  caused  by  an  arbitrary  array  of  surface 
point  charges,  for  carriers  in  a  semiconductor.  We  expressed 
w  in  terms  of  the  statistical  structure  factor  S(Q)  of 
the  surface  charge  array,  using  the  surface  response  function6 
which  gives  the  scattering  potential  of  the  array  as  a  result 
of  screening,,  dielectric,  and  image  effects.  We  show^on 
the  basis  of  statistical  thermodynamic  ensemble  theory,  that 
S(Q)  also  determines  the  chemical  potential  g  of  the  absorbed 
species.  This  dual  role  of  S((£)  enables  us  to  derive  a 
new  and  general  relation  between  surface  scattering  rates 
and  the  adsorption  isotherm.  Transport  measurements  may 
thereby  be  used  to  clarify  chemisorption  processes  and  vice 
versa. 
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I  .  STATISTICAL  STRUCTURE  FACTOR  OF  CHEMISORBED  ATOMS 
Chemisorbed  atoms  make  their  presence  felt  in  carrier 
scattering  through  their  charge,  of  course,  but  also  through 
their  arrangement,  since  that  determines  electron  wave 
interference  effects.  For  scattering  purposes  the  adion 
arrangement  is  best  expressed  by  the  statistical  structure 
factor  S(g),  which  we  will  develop  below.  Having  done  so, 
we  will  then  show  that  S(Q)  can  be  expressed  in  thermodynamic 
terms,  e.g.  the  chemical  potential  of  the  chemisorbed  species 
at  least  m  the  case  where  the  chemisorption  has  reached 
equilibrium. 

An  array  of  point  charges  e  on  a  surface  corresponds  to 
a  surface  charge  density  ft(R): 

-  e  l  «CR  -  5.1 )  (1) 

which  has  the  spatial  average 
Si  =  lim  L 

L-voo 

(flLCR)  is  obtained  by  truncating  ft(R)  to  zero  outside  a 
macroscopic  area  L2.)  We  use  the  autocorrelation  functions 
of  ft  and  6ft  =  ft  -  ft; 


f  dR  aLCR)  (2) 


*.  i.u 


In  scattering  theory  it  is  convenient  to  work  with  the 
Fourier  transforms 

0(Q)  =  (2tt) -1  j dR  fl(R)  e~  i2*  « 


60  (Q) 


=  (  2tt)  " 1  j  dR  60  (R)  e_1S‘S 


=  0(Q)  -  2tt06  (Q)  etc. 


Using  the  truncated  functions  0.  (R)  ,  etc.  one  can 


form  the  spectral  density 


=  lim  L'2|fiT  (g)  ,2 


(4a) 


(4b) 


which  has  the  transform  (Wiener-Khintchine  formula7) 


^n(Q)  =  (4n  2) " 1  J  dAR  n(R)  o(R+AR)  e  ‘  ^ .  (6) 


Likewise,  the  spectral  density  of  60  is 

$6a(S)  *  C2it.!  2  /  dAjR  (D(R)  !!(R+4R)  -  S') 

=  &  (g)  -  S!26  (^)  .  ( 

NM 

We  confine  our  attention  to  chemisorption  charges  e 
which  lie  at  lattice  points  A  which  have  some  definite 
relation  to  the  3-dimensional  crystal  lattice.  Then  Eq.  (1) 
can  be  written 


0(R)  -  e  En(A)6(R  -  A);  (n(A)  =  0,  1  only,)  (8) 

zK 

AU 

where  n(A)  is  the  occupation  number  for  surface  lattice  site  A, 


Of  course, 


0  =  e  n/aQ  , 


-  ( 
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where  n  is  the  mean  occupation  number,  and  aQ  is  the  unit 
cell  area  of  the  surface  lattice.  Using  Eqs.  (6)  and  (8) 
one  can  write 


(eZ/4r2) 


lim  L 
L-*“ 


(L)  CL) 

E  E  m  (A)  n  (A ' ) 
A  A' 

AA 


eiQ.*  C4-A') 


t 


Cio) 


where  the  sum  goes  over  the  points  A^  and  A'  within  the 
area  L“. 

It  is  sufficient  to  consider  only  those  ensembles  which 
have  spatial  uniformity,  in  the  sense  that  (denoting  ensemble 
averages  {}) 


{n (A) n (A+ AA) 1  =  n  p(AA),  independent  of  A.  CH) 

p ( AA)  is  then  the  conditional  occupation  number  of  site  A+AA 

when  n(A)  =  1,  an  p(0)  s  1.  When  Eq.  (11)  is  ensemble 

averaged,  the  E  si  mmation  becomes  trivial  and  one  gets: 

A 

AM 

{$0(Q) }  =  (ne2/4TT2aQ)  Z  p(AA)  e_19  (12a) 


AA 


=  (ne2/4rr2a  )  Z  (p (AA)  -n) e _12  ^  +  U2  l  $  (g-Q)  , 


AA 


2 


(12b) 


where  the  Q  are  all  the  reciprocal  lattice  vectors 

•v. 

Q-A  =  2t r  integer  or 
and  we  used  the  identity 


zero . 


(12c) 


Z  eiS‘A-  =  ( 4 tt 2 / a  )  £,(g  -  £)  Cl2d) 

AA  5 

Instead  of  the  periodic  function  {<^q(Q)),  it  proves  more 
convenient  in  scattering  theory  to  use  the  non-periodic 
spectral  density  of  512  =  Q  -  ft.  From  Eq.  (7)  one  gets 
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l4,CQ)}  =  (ao  &IU2)  s«)  • 


(13a) 


where  S(Q)  is  the  statistical  structure  factor: 

S  CQ)  =  X  ( (p  (AA)  /  n  -  1))  e_1S'A-  +  (47i2/a  )  Z  6(Q-Q).  (13b) 
AA  0  Q^O  v  " 

The  statistical  structure  factor  S(Q)  shows  "Bragg  peaks" 

A 

Q=g  which  correspond  to  complete  order,  n  =  1.  When  n  f  1, 
there  is  generally  also  a  continuous  "background"  which  is 
determined  by  the  pair  correlation  function  p(AA),  i.e.  the 
degree  of  order. 

Having  expressed  the  adion  arrangement  in  terms  of  the 
statistical  structure  factor,  the  form  most  suitable^  for 
calculation  of  carrier  scattering,  we  want  to  bring  out  a 
different  aspect  of  S (Q) .  We  will  show,  in  fact,  that  S (Q) 
can  be  expressed,  by  the  methods  of  statistical  thermodynamics, 
in  terms  of  thermodynamic  properties  chemisorbing  species, 
e.g.  the  chemical  potential  g.  To  do  this  we  have  to  fill 
out  the  description  of  the  system  a  little:  we  let  the  system 
consist  of  chemisorbed  atoms,  of  a  particular  species,  on  a 
crystal  face.  The  system  is  in  thermodynamic  equilibrium 
with  a  reservoir  of  that  species  with  temperature  T  and 
chemical  potential  g.  The  grand  canonical  ensemble  is  then 

O 

appropriate  and  gives  the  exact  result 

nS (0)  =  Z  (p (AA)  -  n)  =  ({N2}-{N}2)/  N  (14) 

AA 

where  N  is  the  fluctuating  total  number  of  atoms  chemisorbed 

on  a  macroscopic  ‘rystal  surface.  The  virtue  of  Eq.  (14)  is 

that  the  R.H.S.  is  an  equilibrium  extensive  parameter 

q 

fluctuation  which  can  be  expressed  exactly  in  terms  of  the 
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•A  (^) 


corresponding  free  energy  properties: 


{N^}  -  { N } ^  -  kT  3{N3-/9g)^,  .  g  =  chemical  potential  (15) 

(The  unit  cell  dimension  aQ  is  considered  as  a  fixed 
extensive  parameter.)  Thus  one  obtains  the  basic  thermo¬ 
dynamic  property  of  the  statistical  structure  factor. 

(S  (0)  }  =  (kT/n2) 3n/3g)T  A  (16) 

’  o 

Some  insight  into  Lq.  (16)  is  provided  by  the  simplest 
adsorption  model,  in  which  interaction  between  adsorbed 
molecules  is  ignored.  The  system  energy  and  grand  partition 
function  are  then 


E  =  ?:  n (A) c  (A)  =  eQN  ;  n(A)  =  0,  1  (17a) 

Q  =  T  exp  ((eQ  -  g)  N/kT)  (17b) 


Because  of  the  formal  similarity  of  this  model  to  that  of  the 
non-interacting  Fermi  gas  one  can  immediately  write  down 

n  =  [1  *  exp((eo-g)/kT)2]'1.  (18) 

This  can  be  seen"*"  to  be  the  familiar  Langmuir  isotherm'*''*'  if 
we  assume  that  the  gas,  with  which  the  adsorbed  species  is 
in  equilibrium,  is  at  moderate  temperature  T  and  pressure  p: 
Then  the  ideal  gas  relation 


e"g/kT  =  (2TrmkT/h2)3/2  (kT/p) 


can  be  used.  This  isotherm  satisfies  the  relation 


kT  Sn 

ry  =  1/n  -  1  (Random  case) 

n  3  g  T 


(19) 


(20) 
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the  other  hand,  the  lack  of  interaction  between  adatoms 
means  we  can  assume  zero  positional  autocorrelation: 
p(AA)  =  n  for  AA  f  0.  Then  Eq.  (13b)  tells  us  that 

SCO)  =  l/n  -  1  (Random  case)  (21) 

Thus  Eq.  (16)  is  satisfied  for  this  simplest  model, 
independent  of  the  general  thermodynamic  fluctuation 
argument . 

It  is  important  to  note  that  Eq.  (20)  breaks  down 
even  for  the  simplest  model  involving  Coulomb  or  other 
interactions.  The  system  energy  can  be  written  as 

E  =  E  n(A)(eo  +  1/2  e0(A));  n(A)  =  0,1  (22) 

A 

where  <f>  (A)  is  the  potential  at  A  due  to  the  other  adions, 
screening ,  and  polarization,  and  hence  depends  on  the  entire 
set  of  occupation  numbers  n(A).  We  make  the  drastic  approx¬ 
imation  of  replacing  0(A)  by  <f>  ,  the  surface  potential  in 

P  ^ 

the  planar  semiconductor  band-bending  model  of  Kingston  and 

Neustadter  in  which  0pg  is  determined  entirely  by  n,  T,  and 

material  constants.  The  grand  partition  function  then 

reduces  to  Eq;  (17b),  except  that  g  is  replaced  by  g  -  l/2ed»  (n)  . 

ps 

We  can  therefore  immediately  write  down  the  "band-bending 
isotherm" 

n  =  [1  +  exp ( (e Q  -  g  +  l/2e0ps  (n) ) /kT] _1  (23) 

in  which  the  coverage  n  is  smaller  than  in  the  Langmuir 
isotherm.  Eq.  (20)  is  now  replaced  by 
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(24) 


77  ~i  =  d/n  -  ID  /  Cl  +  1/2  n2  3u  /3n  T) 
n  3g  ps  i 

IT 

where  u^s  =  e4>pS/kT*  We  cannot.,  howevei ,  conclude  anything 
about  p(AA)  from  such  a  sketchy  treatment  except,  of  course, 
that  Eq.  (16)  must  be  satisfied  because  of  the  general 
equilibrium  fluctuation  arguments. 
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III.  SURFACE  SCATTERING-CHEMICAL  POTENTIAL  RELATION 
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Recently  a  calculation  has  been  made  of  the  surface 
scattering  conduction  electrons  by  a  partially  correlated 
array  of  surface  point  charges,  or,  more  precisely,  by  the 
screened  potential  of  this  array.  The  differential  scattering 
probability  w  for  the  process  |K,k  >  ->•  I K+Q,k * >  (where  K  is  the 

■te*  Z  W.  <X*  Z 

2  2  2 

surface  component  of  the  wave  vector,  and  K  +k  =k'  +(K+ 

was  found,  in  the  Born  approximation,  to  have  the  form 

w(y,n)  =  w(y,n)G0  B(exp  -  2kZj)  S(Q  -  K)  (25) 

Here  wGQ  is  the  scattering  rate  for  a  random  array  of  Yukawa 

potentials  on  the  surface,  B  is  the  surface  space  charge  form 

factor  which  corrects  for  the  actual  form  of  the  space  charge 

potential,  and  (exp  -  2kZj>  is  a  correction  factor  for  the 

dielectric  image  potential.  We  emphasize  that  all  the 

information  about  the  actual  distribution  of  the  scattering 

centers  over  the  surface  is  contained  in  S,  and  the  rest 

of  the  expression  need  be  calculated  only  once  and  applies 

to  any  surface  ulstribution . 

* 

Frequently  in  the  chemisorption  of  ionizing  atoms  on 
semiconductors  the  autocorrelation  p(AA_)  has  a  short  range 
compared  to  the  electron  wavelength: 

p(AAj  =  n  for  all  |  A  A  |  k0  >  1 .  (26) 

In  that  case  we  immediately  get  that 

S(Q  -  K)  =  S(0)  (27) 

<V*A.  -WA 


1  n 


—  (4) 


O') 


l0r  a11  values  of  H  appearing  in  Eq.  (25).  One  then  has, 
using  Eq.  (16), 


w(u,n)  -  w(^n)random  (kT/n(l-n) )  (9n/3g)T  (28) 

which  is  the  fundamental  relation  between  surface  scattering 
rate  and  chemical  potential  of  the  scatterer. 

We  may  now  make  a  direct  and  useful  connection  with  surface 
transport  theory,  e.g.  with  the  semiconductor  field  effect  or 
the  s^ze  effect  of  thin  films.  Here  the  surface  scattering  rates 
enter  the  theory  through  the  Fuchs  reflectivity  p,  according 
to  the  recipe1^ 


p(u)  =  1 


f°  f2 

I  dAr)  (~u/u_  )w(y_  ,n  )  j  1  -  cos 


(29) 


Since  S(0)  is  independent  of  angle,  we  immediately  obtain  an 
analog  of  Eq „  (28)  for  the  Fuchs  reflectivity: 

p(v>  =  P(ll)random  (W/n(l-n»  (9H/8g)T  (30) 

^ow  ^random  can  calcudated^  once  and  for  all  for  a  given 
* 

surface:  it  does  not  depend  upon  the  actual  distribution  of 
the  adions,  but  only  on  their  total  number.  For  any  actual 
chemisorption  equilibrium,  departures  from  the  Langmuir 
isotherm  will  be  produced  by  adion-adion  interactions, 
as  noted  in  the  discussion  of  Eqs .  (20)  and  (24).  These 
departures  can  be  measured,  e.g.,  by  simultaneous  measurements 
of  pressure  and  Hall  coefficient,  the  latter  providing  a  means12* 
of  inferring  the  coverage  n.  The  Fuchs  reflectivity  can  be 
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n 

e*s 


Inferred  from  surface  mobility  measurements.  To  test  Eq.  (30) 
it  would  thus  seem  necessary  to  make  simultaneous  measurements 
of  pressure,  coverage,  and  galvanomagnetic  response  on  a 
single-crystal  semiconductor  face  which  is  otherwide  clean 
to  ultrahigh  vacuum  standards.  Such  measurements,  involving 
various  gases  on  epitaxial  IV-VI  semiconductor  faces,  are  being 
carried  out  at  the  Naval  Ordnance  Laboratory  by  R.  N.  Lee.^ 
Although  the  present  treatment  is  restricted  to  monolayer 
adsorption  with  only  short  range  order,  it  is  hoped  that 
Eq.  (30)  will  provide  a  means  of  applying  observations  of 
transport  inside  a  semiconductor  to  measurements  of  chemisorption 
on  the  outside,  and  vice  versa. 


:.VS  (2>) 
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APPENDIX  A 


The  analog  of  Eq.  (16)  In  the  theory  of  liquids  is  derived^ 
in  ensemble  theory  by  approximating  the  sum  over  states  by  a 
classical  integral  over  3N  particle  coordinates  and  momenta, 
corrected  by  a  factor  (n!h^N)  1  to  account  for  indistinguish- 
ability  and  the  uncertainty  principle.  We  show  below  that 
Eq.  ( 1 6 )  is  exact  for  a  lattice  gas,  by  using  the  occupation 
number  description,  which  is  a  proper  generic  phase  description. 

The  result  is  slightly  different  from  the  classical  continuum 
case . 

From  the  general  weight  functions  PN(n)  and  P(N,n),  for  the 

canonical  and  the  grand  canonical  ensemble,  resp.,  we  want  to 

get  the  two-site  joint  occupation  probabilities  Pn(2)(n(An),  n(A  )) 

( 2 )  N  ~-l  ’  -2 

and  P  (n(A^),  (subscript  N  indicates  the  constraint 

I  n(A)  =  N  =  fixed,  and  5  stands  for  a  "value"  of  the  set  of 
£- 

all  occupation  numbers  n(A)).  One  has 

PN<2)  n(A  ))  =  I '  ' 1  PN(g)  ,  (Al) 

n 

where  the  summation  I'"  is  subject  to  three  constraints: 


n^l^  fixed,  n(A„)  =  fixed,  and  E  n(A)  =  N  =  fixed,  Th^r. 


A 


summing  over  all  values  of  n(A^)  and  n(Ap),  one  gets 


n(A^)  n(A2)  n(*l)a(A2)PN2)(n(A1),n(A2))  =  Z '  n(A1)n(A^,)  I-  (g),  (A2) 

where  E '  is  subject  only  to  the  constraint  Z  n(A)  =  H  =  fixed 
2  A 


Since  n(A)  =  0  or  1  only,  Eq.  (A2)  reduces  t< 
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-  (3) 


PN  (n(.£l}  =  1»  n(A2)  =  !)  =  s  n(A2)  n(A2)  PN(n)  . 


(A3) 


Next,  we  sum  over  all  pairs  of  sites  A^  ^  A^,  getting 


E  E  PN(2)(n(A1)  =  1,  n(A2)  =  1)  =  E  E  E  n(A1)n(A2)  PN(g) 


A,  *  A0 


5  hi  *  *2 


=  N(N-l)  E  PN(g)  -  N(N-l) 
n 


( A-4 ) 


This  procedure  can  be  repeated  for  the  grand  canonical 
ensemble,  using  the  easily  obtained  identity. 

P(2)(n(A1)  =  l,n(A2)  =  D  =  E  P(N)  P^f2)(n(A1)  =  1,  n(A2)  =  1), 


(A5) 


where  P(N)  is  the  probability,  in  the  grand  canonical  ensemble, 
that  there  are  exactly  N  particles  in  the  system.  Summing,  again, 
over  all  pairs  of  sites  A±  /  A2,  one  obtains  from  Eq.  ( A 4 )  that 

E  E  P(2)(n(A, )  =  1,  n(A3)  =  1)  =  E  P(N)  N(N-l)  =  {N2}  -  {N}. 


A-.  ^  A„ 

w-i. 


But,  making  the  identification 


(A6) 


P ^ 2 ^ ( n ( A 1 )  =  1,  n(A2)  =  1)  =  {n( A1)n(A2)}  =  n  p(A2  -  A2),  (A7) 

where  in  the  last  step  we  again  assumed  spatial  uniformity  in 


the  ensemble,  Eq.  (A6)  gives 


{NJ  E  p  ( AA )  =  Or)  -  { N}  , 

AA 

because  {N)/n  is  the  total  number  of  sites.  We  thus  get 


(A8) 
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(A9a) 


{N}  Z  p  ( AA )  =  {N2} 
AA 


or  {N}  Z  (p (AA)  -  n)  =  { N2 }  -  {N}2 
AA 

■u. 

or  {N}  Z  (p ( AA)  -  n)  =  { N2 }  -  {N}2  (1-n). 
AA^O 


In  the  classical  treatment  (1-n)  is  replaced  by  unity 


(A9b) 

(A9c) 
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